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^ . Abstract 

It is noted that one is now in possession of a set of facts, which may be viewed as the 
matching pieces of a puzzle; in that all of them can be resolved by just one idea — that 
is grand unification. These include (i) the observed family-structure, (ii) quantization of 
electric charge, (hi) the meeting of the three gauge couplings, (iv) neutrino oscillations [in 
particular the value Am 2 (^ — i> T ), suggested by SuperK], (v) the intricate pattern of the 
masses and mixings of the fermions, including the smallness of V c b and the largeness of 0° s ^ Vt , 
and (vi) the need for B-L as a generator to implement baryogenesis (via leptogenesis) . All 
these pieces fit beautifully together within a single puzzle board framed by super symmetric 
unification, based on either SO(10) or a string-unified G(224)-symmetry. The two notable 
pieces of the puzzle still missing, however, are proton decay and supersymmetry. 

A concrete proposal is presented within a predictive SO(10)/G(224)-framework that suc- 
cessfully describes the masses and mixings of all fermions, including the neutrinos — with 
eight predictions, all in agreement with observation. Within this framework, a systematic 
study of proton decay is carried out, which (a) pays special attention to its dependence on 
the fermion masses, and (b) limits the threshold corrections so as to preserve natural cou- 
pling unification. The study updates prior work by Babu, Pati and Wilczek, in the context 
of both MSSM and its (interesting) variant, the so-called ESSM, by allowing for improved 
values of the matrix elements and of the short- and long-distance renormalization effects. It 
shows that a conservative upper limit on the proton lifetime is about (1/3 - 2)xl0 34 years, 
with VK + being the dominant decay mode, and quite possibly fi + K° and e + 7r° being promi- 
nent. This in turn strongly suggests that an improvement in the current sensitivity by a 
factor of five to ten (compared to SuperK) ought to reveal proton decay. Otherwise some 
promising and remarkably successful ideas on unification would suffer a major setback. For 
comparison, some alternatives to the conventional approach to unification pursued here are 
mentioned at the end. 
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1 Introduction 



The standard model of particle physics, based on the gauge symmetry SU(2) L xU(l)y x SU(3)c 
]l], is in excellent agreement with observations, at least up to energies of order 100 GeV. 
Its success in turn constitutes a triumph of quantum field theory, especially of the notions of 
gauge invariance, spontaneous symmetry breaking, and renormalizability. The next step in the 
unification-ladder is associated with the concept of "grand unification", which proposes a unity 
of quarks and leptons, and simultaneously of their three basic forces: weak, electromagnetic and 
strong || [|, ||. This concept was introduced on purely aesthetic grounds, in fact before any 
of the empirical successes of the standard model was in place. It was realized in 1972 that the 
standard model judged on aesthetic merits has some major shortcomings || f|]. For example, 
it puts members of a family into five scattered multiplets, assigning rather peculiar hypercharge 
quantum numbers to each of them, without however providing a compelling reason for doing 
so. It also does not provide a fundamental reason for the quantization of electric charge, and it 
does not explain why the electron and proton possess exactly equal but opposite charges. Nor 
does it explain the co-existence of quarks and leptons, and that of the three gauge forces — weak, 
electromagnetic and strong — with their differing strengths. 

The idea of grand unification was postulated precisely to remove these shortcomings. It in- 
troduces the notion that quarks and leptons are members of one family, linked together by a 
symmetry group G, and that the weak, electromagnetic and strong interactions are aspects of 
one force, generated by gauging this symmetry G. The group G of course inevitably contains 
the standard model symmetry G(213) = SU(2)l x U(l)y x SU(3)c as a subgroup. Within this 
picture, the observed differences between quarks and leptons and those between the three gauge 
forces are assumed to be low-energy phenomena that arise through a spontaneous breaking of 
the unification symmetry G to the standard model symmetry G(213), at a very high energy scale 
M 3> 1 TeV. As a prediction of the hypothesis, such differences must then disappear and the true 
unity of quarks and leptons and of the three gauge forces should manifest at energies exceeding 
the scale M. 

The second and perhaps the most dramatic prediction of grand unification is proton decay. 
This important process, which would provide the window to view physics at truly short distances 
(< 10~ 30 cm), is yet to be seen. Nevertheless, as I will stress in this talk, there has appeared over 
the years an impressive set of facts, favoring the hypothesis of grand unification which in turn 
suggest that the discovery of proton decay should be imminent. These include: 

(a) The observed family structure: The five scattered multiplets of the standard model, 
belonging to a family, neatly become parts of a whole (a single multiplet), with their weak hyper- 
charges predicted by grand unification, precisely as observed. It is hard to believe that this is just 
an accident. Realization of this feature calls for an extension of the standard model symmetry 
G(213) = SU(2) L xU(l) y xSU(3) c minimally to the symmetry group G(224) = SU(2) L xSU(2) /? x 
SU(4) C ' ||, which can be extended further into the simple group SO(10) 0, but not SU(5) ||]. The 
G(224) symmetry in turn introduces some additional attractive features (see Section ^), including 
especially the right-handed (RH) neutrinos (vrs) accompanying the left-handed ones (Vl's), and 
B L as a local symmetry. As we will see, both of these features now seem to be needed, on 
empirical grounds, to understand neutrino masses and to implement baryogenesis. 

(b) Quantization of electric charge and the fact that Q e iectron = — Q proton : Grand 
Unification provides compelling reasons for both of these facts. 

(c) Meeting of the gauge couplings: Such a meeting is found to occur at a scale Mx ~ 
2 x 10 16 GeV, when the three gauge couplings are extrapolated from their values measured at 
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LEP to higher energies, in the context of supersymmetry [[7] . This dramatic phenomenon provides 
a strong support in favor of the ideas of both grand unification and supersymmetry || . Both of 
these features in turn may well emerge from a string theory || or M-theory |nj (see discussion 
in Section ||). 

(d) A.m 2 (i/^ — i/ T ) ~ (l/20eV) 2 : The recent discovery of atmospheric neutrino-oscillation 
at SuperKamiokande |Tl| suggests a value Am 2 (i/ M z/ T ) ~ (1/20 eV) 2 . It has been argued (see e.g. 
Ref. |12[) that precisely such a magnitude of Am 2 (z/ M z/ T ) can be understood very simply by 



utilizing the SU(4)-color relation m(z/ r )Dirac ~ ^to P and the SUSY unification scale Mx, noted 
above (See Section [|). 

(e) Some intriguing features of fermion masses and mixings: These include: (i) the 
"observed" near equality of the masses of the b-quark and the r-lepton at the unification-scale 
(i.e. m° ~ mf) and (ii) the observed largeness of the z/^-z/ T oscillation angle (sin 2 20™° Vt > 0.92) 
IPH , together with the smallness of the corresponding quark mixing parameter V c fc(~ 0.04) [[0|]. 
As shown in recent work by Babu, Wilczek and me ||14|| , it turns out that these features and 



more can be understood remarkably well (see discussion in Section |5|) within an economical and 
predictive SO(10)-framework based on a minimal Higgs system. The success of this framework 
is in large part due simply to the group-structure of SO(10). For most purposes, that of G(224) 
suffices. 

(f ) Baryogenesis: To implement baryogenesis [HJ successfully, in the presence of electroweak 



sphaleron effects which wipe out any baryon excess generated at high temperatures in the 
(B-L)-conserving mode, it has become apparent that one would need B-L as a generator of the 
underlying symmetry in four dimensions, whose spontaneous violation at high temperatures would 
yield, for example, lepton asymmetry (leptogenesis). The latter in turn is converted to baryon- 
excess at lower temperatures by electroweak sphalerons. This mechanism, it turns out, yields even 
quantitatively the right magnitude for baryon excess [T7|]. The need for B-L, which is a generator 
of SU(4)-color, again points to the need for G(224) or SO(10) as an effective symmetry near the 
unification-scale Mx- 

The success of each of these six features (a)-(f) seems to be non-trivial. Together they 
make a strong case for both the conventional ideas on supersymmetric grand unification and 
simultaneously for the G(224)/SO(10)-route to such unification, as being relevant to nature at 
short distances < (10 16 GeV) -1 , in four dimensions.^ However, despite these successes, as long 
as proton decay remains undiscovered, the hallmark of grand unification — that is quark-lepton 
trans ■formability — would remain unrevealed. 

The relevant questions in this regard then are: What is the predicted range for the lifetime 
of the proton — in particular an upper limit — within the empirically favored route to unification 
mentioned above? What are the expected dominant decay modes within this route? Are these 
predictions compatible with current lower limits on proton lifetime mentioned above, and if so, 
can they still be tested at the existing or possible near-future detectors for proton decay? 

Fortunately, we are in a much better position to answer these questions now, compared to a 
few years ago, because meanwhile we have learnt more about the nature of grand unification, and 
also there have been improved evaluations of the relevant matrix elements and short and long- 
distance renormalization effects. As noted above (see also Section g and Section the neutrino 
masses and the meeting of the gauge couplings together seem to select out the supersymmetric 
G(224)/SO(10)-route to higher unification. The main purpose of my talk here will therefore be 



1 For comparison, some alternative attempts, including those based on the ideas of (a) large extra dimensions, 
and (b) unification occurring only in higher dimensions, are mentioned briefly in Section 6 G. 
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to address the questions raised above, in the context of this route. For the sake of comparison, 
however, I will state the corresponding results for the case of supersymmetric SU(5) as well. 

My discussion will be based on a recent study of proton decay by Babu, Wilczek and me [fLifl , 
an update presented in the Erice talk flIB| , and a subsequent update of the same as presented 
here. Relative to other analyses, this study has four distinctive features: 

(i) It systematically takes into account the link that exists between proton decay and the 
masses and mixings of all fermions, including the neutrinos. 

(ii) In particular, in addition to the contributions from the so-called "standard" d = 5 op- 
erators |n§ (see Section H), it includes those from a new set of d = 5 operators, related to the 



Majorana masses of the RH neutrinos [[2(J. These latter are found to be generally as important 
as the standard ones. 

(iii) As discussed in the Appendix, the work also restricts GUT-scale threshold corrections, 
so as to preserve naturally coupling unification, in accord with the observed values of the three 
gauge couplings. 

(iv) Finally, the present update incorporates recently improved values of the matrix elements, 
and the short and long-distance renormalization effects, which significantly enhance proton decay 
rate. 

Each of these features turn out to be crucial to gaining a reliable insight into the nature of 
proton decay. Our study shows that the inverse decay rate for the I>K + -mode, which is dominant, 
is less than about 1.2 x 10 31 years for the case of MSSM embedded in minimal SUSY SU(5), and 
that it is less than about 10 33 years for the case of MSSM embedded in SO (10). These upper 
bounds are obtained by making generous allowance for uncertainties in the matrix element and 
the SUSY-spectrum. Typically, the lifetime should of course be less than these bounds. 

Proton decay is studied also for the case of the extended supersymmetric standard model 
(ESSM), that has been proposed a few years ago |2TJ on several grounds, based on the issues 
of (a) an understanding of the inter-family mass- hierarchy, (b) removing the mismatch between 
MSSM and string- unification scales, and (c) dilaton-stabilization (see Section ^ and the appendix). 
This case adds an extra pair of vector-like families at the TeV-scale, transforming as 16 + 16 of 
SO(10), to the MSSM spectrum. While the case of ESSM is fully compatible with both neutrino- 
counting at LEP and precision electroweak tests, it can of course be tested directly at the LHC 
through a search for the vectorlike fermions. Our study shows that, with the inclusion of only 
the "standard" d = 5 operators (defined in Section [|), ESSM, embedded in SO(10), can quite 
plausibly lead to proton lifetimes in the range of 10 33 — 10 34 years, for nearly central values of 
the parameters pertaining to the SUSY-spectrum and the matrix element. Allowing for a wide 
variation of the parameters, owing to the contributions from both the standard and the neutrino 
mass-related d = 5 operators (discussed in Section ||), proton lifetime still gets bounded above by 
about 2 x 10 34 years, for the case of ESSM, embedded in SO(10) or a string-unified G(224). 

For either MSSM or ESSM, embedded in G(224) or SO(10), due to contributions from the 
new operators, the n + K°-mode is found to be prominent, with a branching ratio typically in the 
range of 10-50%. By contrast, minimal SUSY SU(5), for which the new operators are absent, 
would lead to branching ratios < 10~ 3 for this mode. It is stressed that the e + 7r°-mode induced 
by gauge boson-exchange, in either SUSY SU(5) or SUSY SO(10), could have an inverse decay 
rate as short as about (1 — 2) x 10 34 years. 

Thus our study of proton decay, correlated with fermion masses, strongly suggests that discov- 
ery of proton decay should be imminent. Allowing for the possibility that the proton lifetime may 
well be closer to the upper bound stated above, a next-generation detector providing a net gain 
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in sensitivity in proton decay-searches by a factor of 5-10, compared to SuperK, would certainly 
be needed not just to produce proton-decay events, but also to clearly distinguish them from 
the background. It would of course also be essential to study the branching ratios of certain sub- 
dominant but crucial decay modes, such as the /j, + K° and e + 7r°. The importance of such improved 
sensitivity, in the light of the successes of supersymmetric grand unification, is emphasized at the 
end. 



2 Advantages of the Symmetry G(224) as a Step to Higher Unifica- 
tion 

As mentioned in the introduction, the hypothesis of grand unification was introduced to remove 
some of the conceptual shortcomings of the standard model (SM). To illustrate the advantages of 
an early suggestion in this regard, consider the five standard model multiplets belonging to the 
electron-family as shown: 
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Here the superscripts denote the respective weak hypercharges Y w (where Q em = I 3 l + Y w /2) 
and the subscripts L and R denote the chiralities of the respective fields. If one asks: how one can 
put these five multiplets into just one multiplet, the answer turns out to be simple and unique. 
As mentioned in the introduction, the minimal extension of the SM symmetry G(213) needed, to 
achieve this goal, is given by the gauge symmetry ||: 

G(224) = SU(2) L x SU(2) fl x SU(4) C . (2) 

Subject to left-right discrete symmetry [L <-> R), which is natural to G(224), all members of the 
electron family become parts of a single left-right self-conjugate multiplet, consisting of: 



F L,R 



U r Uy U b V e 

d r d v db e~ 



(3) 



The multiplets F[ and are left-right conjugates of each other and transform respectively as 
(2,1,4) and (1,2,4) of G(224); likewise for the muon and the tau families. Note that the symmetries 
SU(2) L and SU(2) R are just like the familiar isospin symmetry, except that they operate on quarks 
and well as leptons, and distinguish between left and right chiralities. The left weak-isospin SU(2)i 
treats each column of F[ as a doublet; likewise SU(2)# for Fr. The symmetry SU(4)-color treats 
each row of F[ and F R as a quartet; thus lepton number is treated as the fourth color. Note also 
that postulating either SU(4)-color or SU(2)# forces one to introduce a right-handed neutrino (ur) 
for each family as a singlet of the SM symmetry. This requires that there be sixteen two-component 
fermions in each family, as opposed to fifteen for the SM. The symmetry G(224) introduces an 
elegant charge formula: 

Qem = hh + hR + g ( 4 ) 

expressed in terms of familiar quantum numbers I^l-, I3R and -B-L, which applies to all forms of 
matter (including quarks and leptons of all six flavors, gauge and Higgs bosons). Note that the 
weak hypercharge given by Yw/2 = I^r + 3 jT L is now completely determined for all members of 
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the family. The values of IV thus obtained precisely match the assignments shown in Eq. (]T|). 
Quite clearly, the charges I 3 r and B L, being generators respectively of SU(2) L , SU(2)# and 
SU(4) C , are quantized; so also then is the electric charge Q em . 

In brief, the symmetry G(224) brings some attractive features to particle physics. These 
include: 

(i) Unification of all 16 members of a family within one left-right self-conjugate multiplet; 

(ii) Quantization of electric charge, with a reason for the fact that ^electron = — Qproton 

(iii) Quark-lepton unification (through SU(4) color); 

(iv) Conservation of parity at a fundamental level || |22| ; 

(v) Right-handed neutrinos {v' R s) as a compelling feature; and 
fvi) B-L local symmetry. 

As mentioned in the introduction, the two distinguishing features of G(224) — i.e. the existence of 
the RH neutrinos and B-L as a local symmetry — now seem to be needed on empirical grounds. 
Furthermore, SU(4)-color provides simple relations between the masses and mixings of quarks and 
leptons, while SU(2) L x SU(2) fi relates the mass-matrices in the up and down sectors. As we will 
see in Sections [| and ^, these relations are in good accord with observations. 

Believing in a complete unification, one is led to view the G(224) symmetry as part of a bigger 
symmetry, which itself may have its origin in an underlying theory, such as string theory. In this 
context, one may ask: Could the effective symmetry below the string scale in four dimensions 
(see Section ||) be as small as just the SM symmetry G(213), even though the latter may have its 
origin in a bigger symmetry, which lives only in higher dimensions? I will argue in Section [| that 
the data on neutrino masses and the need for baryogenesis provide an answer to the contrary, 
suggesting that it is the effective symmetry in four dimensions, below the string scale, which must 
minimally contain either G {22 4) or a close relative G(214) = SU(2) i xI 3J jxSU(4) c '. 

One may also ask: does the effective four dimensional symmetry have to be any bigger than 
G(224) near the string scale? In preparation for an answer to this question, let us recall that 
the smallest simple group that contains the SM symmetry G(213) is SU(5) It has the virtue 
of demonstrating how the main ideas of grand unification, including unification of the gauge 
couplings, can be realized. However, SU(5) does not contain G(224) as a subgroup. As such, it 
does not possess some of the advantages listed above. In particular, it does not contain the RH 
neutrinos as a compelling feature, and B-L as a local symmetry. Furthermore, it splits members 
of a family (not including vr) into two multiplets: 5 + 10. 

By contrast, the symmetry SO(10) has the merit, relative to SU(5), that it contains G(224) as 
a subgroup, and thereby retains all the advantages of G(224) listed above. (As a historical note, 
it is worth mentioning that these advantages had been motivated on aesthetic grounds through 
the symmetry G(224) ||, and all the ideas of higher unification were in place || ||, ||, before 
it was noted that G(224) [isomorphic to SO (4) x SO (6)] embeds nicely into SO (10) @). Now, 
SO(10) even preserves the 16-plet family-structure of G(224) without a need for any extension. 
By contrast, if one extends G(224) to the still higher symmetry E 6 ||23|| , the advantages (i)-(vi) 
are retained, but in this case, one must extend the family-structure from a 16 to a 27-plet, by 
postulating additional fermions. In this sense, there seems to be some advantage in having the 
effective symmetry below the string scale to be minimally G(224) [or G(214)] and maximally no 
more than SO(10). I will compare the relative advantage of having either a string-derived G(224) 
or a string-SO(lO), in the next section. First, I discuss the implications of the data on coupling 
unification. 
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3 The Need for Supersymmetry: MSSM versus String Unifications 



It has been known for some time that the precision measurements of the standard model coupling 
constants (in particular sin 2 6\y) at LEP put severe constraints on the idea of grand unification. 
Owing to these constraints, the non-supersymmetric minimal SU(5), and for similar reasons, the 



one-step breaking minimal non-supersymmetric SO(10)-model as well, are now excluded [23]. But 
the situation changes radically if one assumes that the standard model is replaced by the minimal 
supersymmetric standard model (MSSM), above a threshold of about 1 TeV. In this case, the 
three gauge couplings are found to meet 0, to a very good approximation, barring a few percent 
discrepancy which can be attributed to threshold corrections (see Appendix). Their scale of 
meeting is given by 

M x ~2x 10 16 GeV (MSSM or SUSY SU(5)) . (5) 

This dramatic meeting of the three gauge couplings, or equivalently the agreement of the 
MSSM-based prediction of sin 9w( m z)tb — 0.2315 ± 0.003 f25| with the observed value of 
sin 2 Owijnz) = 0.23124 ± 0.00017 IT3"] , provides a strong support for the ideas of both grand 
unification and supersymmetry, as being relevant to physics at short distances < (10 16 GeV) -1 . 

In addition to being needed for achieving coupling unification there is of course an independent 
motivation for low-energy supersymmetry — i.e. for the existence of SUSY partners of the standard 
model particles with masses of order 1 TeV. This is because it protects the Higgs boson mass from 
getting large quantum corrections, which would (otherwise) arise from grand unification and 
Planck scale physics. It thereby provides at least a technical resolution of the so-called gauge- 
hierarchy problem. In this sense low-energy supersymmetry seems to be needed for the consistency 
of the hypothesis of grand unification. Supersymmetry is of course also needed for the consistency 
of string theory. Last but not least, as a symmetry linking bosons and fermions, it is simply a 
beautiful idea. And it is fortunate that low-energy supersymmetry can be tested at the LHC, and 
possibly at the Tevatron, and the proposed NLC. 

The most straightforward interpretation of the observed meeting of the three gauge couplings 
and of the scale Mx, is that a supersymmetric grand unification symmetry (often called GUT 
symmetry), like SU(5) or SO(10), breaks spontaneously at M x into the standard model symmetry 
G(213), and that supersymmetry-breaking induces soft masses of order one TeV. 

Even if supersymmetric grand unification may well be a good effective theory below a certain 
scale M > Mx, it ought to have its origin within an underlying theory like the string/M theory. 
Such a theory is needed to unify all the forces of nature including gravity, and to provide a 
good quantum theory of gravity. It is also needed to provide a rationale for the existence of 
flavor symmetries (not available within grand unification), which distinguish between the three 
families and can resolve certain naturalness problems including those associated with inter-family 
mass hierarchy. In the context of string or M-theory, an alternative interpretation of the observed 
meeting of the gauge couplings is however possible. This is because, even if the effective symmetry 
in four dimensions emerging from a higher dimensional string theory is non-simple, like G(224) or 
even G(213), string theory can still ensure familiar unification of the gauge couplings at the string 
scale. In this case, however, one needs to account for the small mismatch between the MSSM 
unification scale Mx (given above), and the string unification scale, given by M st ~ g s t x 5.2 x 10 17 
GeV « 3.6 x 10 17 GeV (Here we have put a st = a GC/T (MSSM) « 0.04) |3§. Possible resolutions 
of this mismatch have been proposed. These include: (i) utilizing the idea of string- duality [p7 



which allows a lowering of M st compared to the value shown above, or alternatively (ii) the idea of 
the so-called "Extended Supersymmetric Standard Model" (ESSM) that assumes the existence of 
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two vector-like families, transforming as (16 + 16) of SO(10), with masses of order one TeV [tH 
addition to the three chiral families. The latter leads to a semi-perturbative unification by raising 
olqut to about 0.25-0.3. Simultaneously, it raises Mx, in two loop, to about (1/2 — 2) x 10 17 GeV. 
(Other mechanisms resolving the mismatch are reviewed in Ref. ||28||). In practice, a combination 
of the two mechanisms mentioned above may well be relevant. Q 

While the mismatch can thus quite plausibly be removed for a non-GUT string-derived sym- 
metry like G(224) or G(213), a GUT symmetry like SU(5) or SO(10) would have an advantage 
in this regard because it would keep the gauge couplings together between M st and Mx (even if 
Mx ~ M st /20), and thus not even encounter the problem of a mismatch between the two scales. 
A supersymmetric four dimensional GUT-solution [like SU(5) or SO(10)], however, has a possible 
disadvantage as well, because it needs certain color triplets to become superheavy by the so-called 
doublet-triplet splitting mechanism (see Section |6| and Appendix), in order to avoid the problem 
of rapid proton decay. However, no such mechanism has emerged yet, in string theory, for the 
GUT-like solutions @. Q 

Non-GUT string solutions, based on symmetries like G(224) or G(2113) for example, have 
a distinct advantage in this regard, in that the dangerous color triplets, which would induce 
rapid proton decay, are often naturally projected out for such solutions 0, [Hp. Furthermore, 
the non-GUT solutions invariably possess new "flavor" gauge symmetries, which distinguish be- 
tween families. These symmetries are immensely helpful in explaining qualitatively the observed 
fermion mass-hierarchy (see e.g. Ref. [|3l|]) and resolving the so-called naturalness problems of 



supersymmetry such as those pertaining to the issues of squark- degeneracy p2f , CP violation pi 



and quantum gravity- induced rapid proton decay |3j . 

Weighing the advantages and possible disadvantages of both, it seems hard at present to make 
a priori a clear choice between a GUT versus a non-GUT string-solution. As expressed elsewhere 



35 1, it therefore seems prudent to keep both options open and pursue their phenomenological 



consequences. Given the advantages of G(224) or SO(10) in the light of the neutrino masses (see 
Sections || and §), I will thus proceed by assuming that either a suitable four dimensional G(224)- 
solution [with the scale Mx being close to M st (see footnote 2)], or a realistic four- dimensional 
SO(10)-solution (with the desired mechanism for doublet-triplet splitting) emerges effectively 
from an underlying string theory, at the "conventional" string-scale M st ~ 10 17 -10 18 GeV, and 
that the G(224)/SO(10) symmetry in turn breaks spontaneously at the conventional GUT-scale 
of M x ~ 2 x 10 16 GeV (or at M x ~ 5 x 10 16 GeV for the case of ESSM, as discussed in footnote 2) 
to the standard model symmetry G(213). The extra dimensions of string/M-theory are assumed 
to be tiny with sizes < Mx ~ 10 -30 cm, so as not to disturb the successes of GUT. In short, 
I assume that essentially the conventional (good old) picture of grand unification, proposed and 
developed sometime ago || ; [|, [|, 0/, holds as a good effective theory above the unification 



2 I have in mind the possibility of string-duality [|27| lowering M s t for the case of semi-perturbative unification 
in ESSM (for which a s t ~ 0.25, and thus, without the use of string-duality, M s t would have been about 10 18 GeV) 
to a value of about (1-2) xlO 17 GeV (say), and semi-perturbative unification ||l| raising the MSSM value of Mx 
to about 5xl0 16 GeVw M st (l/2 to 1/4) (say). In this case, an intermediate symmetry like G(224) emerging at 
M st would be effective only within the short gap between M st and Mx, where it would break into G(213). Despite 
this short gap, one would still have the benefits of SU(4)-color that are needed to understand neutrino masses 
(see Section 4), and to implement baryogenesis via leptogenesis. At the same time, since the gap is so small, 
the couplings of G(224), unified at M s t would remain essentially so at Mx, so as to match with the "observed" 
coupling unification, of the type suggested in Ref. pjj . 

3 Some alternative mechanisms for doublet-triplet splitting, and for suppression of the d = 5 proton decay 
operators have been proposed in the context of higher dimensional theories. These will be mentioned briefly in 
Section 6 G. 
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scale Mx and up to some high scale M < M st , with the added presumption that it may have its 
origin from the string / M-theory . Such a picture seems to be directly motivated on observational 
grounds such as those based on (a) coupling unification (discussed above), (b) neutrino masses 
including the (mass) 2 -difference of the v^-Vt system and the near maximal v il -v T oscillation angle 
(see discussions in the next sections), and (c) the fact that spontaneous violation of B L local 
symmetry at high temperatures, seems to be needed to implement baryogenesis via leptogenesis.Q 
We will see that with the broad assumption mentioned above, an economical and predictive 
framework emerges, which successfully accounts for a host of observed phenomena pertaining 
to the masses and the mixings of all fermions, including neutrinos. It also makes some crucial 
testable predictions for proton decay. I next discuss the implications of the mass of z/ T , or rather 
of Am 2 (i^z/ T ), as revealed by the SuperK data. 



4 Am 2 (z/ M i/ r ): Evidence In Favor of the G(224) Route 

One can obtain an estimate for the mass of v T L in the context of G(224) or SO(10) by using the 
following three steps (see e.g. Ref. Hl2| ): 

(i) Assume that B— L and I3R, contained in a string-derived G(224) or SO(10), break near the 
unification-scale: 

M x ~ 2 x 10 16 GeV, (6) 

through VEVs of Higgs multiplets of the type suggested by string-solutions — i.e. ((1,2,4)//) for 
G(224) or (16//) for SO(10), as opposed to 126// which seems to be unobtainable at least in 
weakly interacting string theory |B9j. In the process, the RH neutrinos which are singlets 
of the standard model, can and generically will acquire superheavy Majorana masses of the type 
Mjj v l fi C~ l u R , by utilizing the VEV of (16//) and effective couplings of the form: 

C M (50(10)) = fa 16, ■ 16, 1Q H ■ 16 H /M + h.c. (7) 

A similar expression holds for G(224). Here i,j = 1,2,3, correspond respectively to e, /1 and 
t families. Such gauge-invariant non-renormalizable couplings might be expected to be induced 
by Planck-scale physics, involving quantum gravity or stringy effects and/or tree-level exchange 
of superheavy states, such as those in the string tower. With (at least the largest among 
them) being of order unity, we would thus expect M to lie between Mpi anc k ~ 2 x 10 18 GeV and 
-^string ~ 4 x 10 17 GeV. Ignoring for the present off-diagonal mixings (for simplicity), one thus 
obtains f\ 

M 3R « /33( ^ h)2 « /aa (2 x 10 14 GeV) p 2 (M Planck /M) (8) 

This is the Majorana mass of the RH tau neutrino. Guided by the value of Mx, we have 
substituted (16//) = (2 x 10 16 GeV) p , where we expect p « 1/2 to 2 (say). 

4 Alternative scenarios such as those based on TeV-scale large extra dimensions (3(| or string-scale being at a 
few TeV f37| , or submillimeter-size even larger extra dimensions with the fundamental scale of quantum gravity 
being a few TeV [Q, though intriguing, do not seem to provide simple explanations of these features: (a), (b) and 
(c). They will be mentioned briefly in Section 6 G. 

5 The effects of neutrino-mixing and of the more legitimate choice of M = M str ing ~ 4 x 10 17 GeV (instead of 
M = A/pi anc k) on the values of ?n(^£) and of M3R are considered in Ref. juj and are reflected in our discussions 
in Section 5. The two effects together end up in yielding essentially the same mass for m{y£) as obtained within 
the simplified picture presented in this section, together with a value for M 3 /j « (5-10) x 10 14 GeV. 
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(ii) Now using SU(4)-color and the Higgs multiplet (2, 2, 1)# of G(224) or equivalently 10h 
of SO(10), one obtains the relation m T (Mx) = m^Mx), which is known to be successful. Thus, 
there is a good reason to believe that the third family gets its masses primarily from the 10# or 
equivalently (2,2, 1)# (see Section 5). In turn, this implies: 

^KiraJ « rn top (M x ) « (100-120)GeV . (9) 

Note that this relationship between the Dirac mass of the tau-neutrino and the top-mass is special 
to SU(4)-color. It does not emerge in SU(5). 

(iii) Given the superheavy Majorana masses of the RH neutrinos as well as the Dirac masses 
as above, the see-saw mechanism j|0| yields naturally light masses for the LH neutrinos. For v T L 



(ignoring flavor- mixing) , one thus obtains, using Eqs.® and (g), 

miyl) » m( ^ irac)2 « [(1/20) eV (1 - 1.44)// 33 p 2 } (M/M Planck ) . (10) 

In the next section, we discuss the masses and mixings of all three neutrinos. As we will 
see, given the hierarchical masses of quarks and charged leptons and the see-saw mechanism, we 
naturally obtain m(y^) ~ (l/10)m(i>£). We are thus led to predict that Am 2 (z/ /J z/ r ) t / l = \vr?{y T L ) 



m 



u L)\th ~ m2 ( zy I = square of the RHS of Eq. (|i~0P . Now SuperK result strongly suggests 
that it is observing ^l'^l (rather than v^-Vx) oscillation, with a 
Aw^^^jois ~ 3x 10~ 3 eV 2 . It seems truly remarkable that the expected magnitude of Ara 2 (^!/ T ), 
given to a very good approximation by the square of the RHS of Eq. (|T0|), is just about what is 
observed at SuperK, if / 33 p 2 (M Planck /M) 1.3 to 1/2. Such a range for / 33 p 2 (M Planck /M) seems 
most plausible and natural (see discussion in Ref. |H§)- Note that the estimate ( PH| ) crucially 
depends upon the supersymmetric unification scale, which provides a value for M 3 #, as well as 
on SU(4)-color that yields ^(^5irac)- ^ e agreement between the expected and the SuperK results 
thus clearly favors supersymmetric unification, and in the string theory context, it suggests that 
the effective symmetry below the string-scale should contain SU(4)-color. Thus, minimally this 
effective symmetry should be either G(214) or G(224), and maximally as big as SO(10), if not Eq. 

By contrast, if SU(5) is regarded as either a fundamental symmetry or as the effective symmetry 
below the string scale, there would be no compelling reason based on symmetry alone, to introduce 
a u R , because it is a singlet of SU(5). Second, even if one did introduce v R by hand, their Dirac 
masses, arising from the coupling h l 5^(5^)^, would be unrelated to the up-flavor masses and 
thus rather arbitrary [contrast with Eq. @]. So also would be the Majorana masses of the fjj's, 
which are SU(5)-invariant, and thus can be even of order string scale . This would give extremely 
small values of m[y r L ) and vniy^) and thus of Am 2 ^), which would be in gross conflict with 
observation. 

Before passing to the next section, it is worth noting that the mass of v T or of Am 2 ^^) 
suggested by SuperK, as well as the observed value of sin 2 6 W (see Section |3|), provide valuable 
insight into the nature of GUT symmetry breaking. They both favor the case of a single-step 
breaking (SSB) of SO(10) or a string-unified G(224) symmetry at a high scale of order Mx, into 
the standard model symmetry G(213), as opposed to that of a multi-step breaking (MSB). The 
latter would correspond, for example, to SO(10) [or G(224)] breaking at a scale M x into G(213), 
which in turn breaks at a scale M 2 -C Mi into G(213). One reason why the case of single-step 
breaking is favored over that of MSB is that the latter can accommodate but not really predict 
sin 2 9w, whereas the former predicts the same successfully. Furthermore, since the Majorana 
mass of v T R arises arises only after B L and J 3 # break, it would be given, for the case of MSB, by 
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M 3R ~ / 33 (Mf/M), where M ~ M st (say). If M 2 < M x ~ 2 x 10 16 GeV, and M > M x , one 
would obtain too low a value (<C 10 14 GeV) for M 3 # [compare with Eq. (8)], and thereby too large 
a value for m(v£), compared to that suggested by SuperK. By contrast, the case of single-step 
breaking (SSB) yields the right magnitude for m{v T ) [see Eq. (10)]. 

Thus the success of the results on m(u T ) and thereby on Am 2 (z/ M z/ r ) discussed above not only 
favors the symmetry SO (10) or G(224) beging effective in J^D at a high scale, but also clearly 
suggests that B-L and I 3 r break near the conventional GUT scale Mx ~ 2 x 10 16 GeV, rather 
than at an intermediate scale <C M x ■ In other words, the observed values of both sin 2 6w and 
Am 2 (^i/ T ) favor only the simplest pattern of symmetry-breaking, for which SO(10) or a string- 
derived G(224) symmetry breaks in one step to the standard model symmetry, rather than in 
multiple steps. It is of course only this simple pattern of symmetry breaking that would be rather 
restrictive as regards its predictions for proton decay (to be discussed in Section ||). 1 nex t discuss 
the problem of understanding the masses and mixings of all fermions. 



5 Understanding Aspects of Fermion Masses and Neutrino Oscilla- 
tions in SO(10) 

Understanding the masses and mixings of all quarks in conjunction with those of the charged 
leptons and neutrinos is a goal worth achieving by itself. It also turns out to be essential for 
the study of proton decay. I therefore present first a partial attempt in this direction, based 



on a quark-lepton unified G(224)/SO(10)-framework, which seems most promising [|14j]. A few 
guidelines would prove to be helpful in this regard. The first of these is motivated by the desire 
for economy [see flnp], and the rest (see below) by the data. In essence, we will be following 
(partly) a bottom-up approach by appealing to the data to provide certain clues as regards the 
pattern of the Yukawa couplings, and simultaneously a top-down approach by appealing to grand 
unification, based on the symmetry G(224)/SO(10), to restrict the couplings by the constraints 
of group theory. The latter helps to interrelate the masses and mixings of quarks with those of 
the charged leptons and the neutrinos. As we will see, it is these interrelationships, which permit 
predictivity, and are found to be remarkably successful. The guidelines which we adopt are as 
follows. 



1) Hierarchy Through Off-diagonal Mixings: Recall earlier attempts |4]J that attribute 
hierarchical masses of the first two families to mass matrices of the form: 

M=(° [}mf\ (11) 

for the (d, s) quarks, and likewise for the (u, c) quarks. Here e ~ 1/10. The hierarchical patterns 
in Eq. (|ll|) can be ensured by imposing a suitable flavor symmetry which distinguishes between 
the two families (that in turn may have its origin in string theory (see e.g. Ref |Jl]). Such a 
pattern has the virtues that (a) it yields a hierarchy that is much larger than the input parameter 
e: {m d /m s ) ~ 6 2 C 6, and (b) it leads to an expression for the Cabibbo angle: 





12) 



which is rather successful. Using ^md/m s rs 0.22 and \jm u jm c ~ 0.06, we see that Eq. (|12|) 
works to within about 25% for any value of the phase 0. Note that the square root formula (like 
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^Jmd/m s ) for the relevant mixing angle arises because of the symmetric form of M in Eq. (|TT|), 
which in turn is ensured if the contributing Higgs is a 10 of SO(10). A generalization of the 
pattern in Eq. ([IT]) would suggest that the first two families (i.e. the e and the /z) receive masses 
primarily through their mixing with the third family (r), with (1,3) and (1,2) elements being 
smaller than the (2, 3); while (2, 3) is smaller than the (3, 3). We will follow this guideline, except 
for the modification noted below. 

2) The Need for an Antisymmetric Component: Although the symmetric hierarchical 
matrix in Eq. ( [TT]) works well for the first two families, a matrix of the same form fails altogether 
to reproduce V c b, for which it yields: 



cb 




(13) 



Given that ^Jm s /rrib « 0.17 and ^m c /m t ~ 0.0.06, we see that Eq. ( |T3| ) would yield V c b varying 
between 0.11 and 0.23, depending upon the phase x- This is too big, compared to the observed 
value of V c b ~ 0.04 ± 0.003, by at least a factor of 3. We interpret this failure as a clue to 
the presence of an antisymmetric component in M, together with symmetrical ones (so that 
rriij 7^ niji), which would modify the relevant mixing angle to \Jrrii/ rrij ^Jrriij/mji, where m; and 
rrij denote the respective eigenvalues. 

3) The Need for a Contribution Proportional to B— L: The success of the relations 
m" ~ m°, and m° ~ ^(^ T )r)irac ( see Section |]), suggests that the members of the third family get 
their masses primarily from the VEV of a SU(4)-color singlet Higgs field that is independent of 
B-L. This is in fact ensured if the Higgs is a 10 of SO(10). However, the empirical observations 
of m Q s ~ ^°/3 and m\ ~ 3m° call for a contribution proportional to B-L as well. Further, 
one can in fact argue that understanding naturally the suppression of V c b (in the quark-sector) 
together with an enhancement of 0° s ^ Ut (in the lepton sector) calls for a contribution that is not only 
proportional to B-L, but also antisymmetric in the family space (this later feature is suggested 
already in item (^)). We show below how both of these requirements can be met in SO(10), even 
for a minimal Higgs system. 

4) Up-Down Asymmetry: Finally, the up and the down-sector mass matrices must not be 
proportional to each other, as otherwise the CKM angles would all vanish. Note that the cubic 
couplings of a single 10// with the fermions in the 16's will not serve the purpose in this regard. 



Following Ref. fl4|| , I now present a simple and predictive mass-matrix, based on SO(10), that 
satisfies all four requirements (|lj), (0), (|3|) and (|). The interesting point is that one can obtain 
such a mass-matrix for the fermions by utilizing only the minimal Higgs system, that is needed 
anyway to break the gauge symmetry SO(10). It consists of the set: 

^minimal = {45//, 16//, 16//, 10//} . (14) 

Of these, the VEV of (45//) ~ M x breaks SO(10) into G(2213), and those of (16//) = (16//) ~ 
Mx break G(2213) to G(213), at the unification-scale Mx- Now G(213) breaks at the electroweak 
scale by the VEV of (10//) to U(l) em x SU(3) C . 

One might have introduced large- dimensional tensorial multiplets of SO(10) like 126// and 
120//, both of which possess cubic level Yukawa couplings with the fermions. In particular, the 
coupling 16jl6j(120//) would give the desired family-antisymmetric as well as (B-L)-dependent 
contribution. We do not however introduce these multiplets in part because there is a general 
argument suggesting that they do not arise at least in weakly interacting heterotic string solutions 
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39| , and in part also because mass-splittings within such large- dimensional multiplets could give 
excessive threshold corrections to a 3 {m z ) (typically exceeding 20%), rendering observed coupling 
unification fortuitous. By contrast, the multiplets in the minimal set (shown above) can arise 
in string solutions. Furthermore, the threshold corrections for the minimal set are found to be 



naturally small, and even to have the right sign, to go with the observed coupling unification 
(see Appendix). 

The question is: can the minimal set of Higgs multiplets [see Eq. (|14|)1 meet all the requirements 
listed above? Now 10# (even several 10's) cannot meet the requirements of antisymmetry and 
(.B-L)-dependence. Furthermore, a single 10 h cannot generate CKM-mixings. This impasse 
disappears, however, as soon as one allows for not only cubic, but also effective non-renormalizable 
quartic couplings of the minimal set of Higgs fields with the fermions. These latter couplings 
could of course well arise through exchanges of superheavy states (e.g. those in the string tower) 
involving renormalizable couplings, and/or through quantum gravity. 

Allowing for such cubic and quartic couplings and adopting the guideline (JIJ) of hierarchical 
Yukawa couplings, as well as that of economy, we are led to suggest the following effective la- 
grangian for generating Dirac masses and mixings of the three families |TJ| (for a related but 



different pattern, involving a non- minimal Higgs system, see Ref. ||43|| ). 

£ Yuk = h 33 16 3 16 3 10h + [ h 23 16 2 16 3 10^ + a 23 16 2 16 3 10^ 45 H /M 
+ g 23 16 2 16 3 16^ 16 H /M] + {a 12 16 1 16 2 10 H 45 // /M 

+ g 12 16i 16 2 16^ 16 H /M} . (15) 

Here, M could plausibly be of order string scale. Note that a mass matrix having essentially 
the form of Eq. (|TTJ) results if the first term /i 33 (10#) is dominant. This ensures « and 
m' m m°(^Dirac)- Following the assumption of progressive hierarchy (equivalently appropriate 
flavor symmetries 0), we presume that h 23 ~ h 33 /10, while h 22 and tin, which are not shown, 
are assumed to be progressively much smaller than h 23 . Since (45#) ~ (16#) ~ M x , while 
M ~ M st ~ 10M X , the terms a 23 (45#)/M and g 23 (16n) /M can quite plausibly be of order 
/?.33/10, if a 23 ~ g 23 ~ h 33 . By the assumption of hierarchy, we presume that a\ 2 <C a 23 , and 

912 < 923 

It is interesting to observe the symmetry properties of the a 2 3 and g 2 3-terms. Although 10 h x 
45h = 10 + 120 + 320, given that (45#) is along B-L, which is used to implement doublet-triplet 
splitting (see Appendix), only 120 in the decomposition contributes to the mass- matrices. This 
contribution is, however, antisymmetric in the family-index and, at the same time, proportional 
to B-L. Thus the a 23 term fulfills the requirements of both antisymmetry and (B-L) -dependence, 
simultaneously [] . With only and a^-terms, however, the up and down quark mass-matrices will 
be proportional to each other, which would yield Vckm = 1- This is remedied by the coupling, 

6 Although no explicit string solution with the hierarchy in all the Yukawa couplings in Eq. (|l5| ) — i.e. in hij, 
ciij and gij — exists as yet, one can postulate flavor symmetries of the type alluded to (e.g. two abelian U(l) 
symmetries), which assign flavor charges not only to the fermion families and the Higgs multiplets, but also to a 
few (postulated) SM singlets that acquire VEVs of order Mx- The flavor symmetry-allowed effective couplings 
such as 16216310^ (S) /M would lead to /123 ~ (S) /M ~ f/fO. One can verify that the full set of hierarchical 
couplings shown in Eq. (|l5|) can in fact arise in the presence of two such U(l) symmetries. String theory (at least) 



offers the scope (as indicated by the solutions of Refs. |3]J and |30|) for providing a rationale for the existence 
of such flavor symmetries, together with that of the SM singlets. For example, there exist solutions with the top 
Yukawa coupling being leading and others being hierarchical (as in Ref. pl|). 

7 The analog of 10# • 45h for the case of G(224) would be xh = (2, 2, 1)h ■ (1, 1, 15)jj. Although in general, the 
coupling of xh to the fermions need not be antisymmetric, for a string-derived G(224), the multiplet (f ,l,15)jj is 
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because, the 16h can have a VEV not only along its SM singlet component (transforming as 
Vr) which is of GUT-scale, but also along its electroweak doublet component — call it 16^ — of the 
electroweak scale. The latter can arise by the the mixing of 16^ with the corresponding doublet 
(call it lOd) in the 10#. The MSSM doublet Hd, which is light, is then a mixture of 10^ and 16^, 
while the orthogonal combination is superheavy (see Appendix). Since (16^) contributes only 
to the down-flavor mass matrices, but not to the up-flavor, the g 23 and g± 2 couplings generate 
non-trivial CKM-mixings. We thus see that the minimal Higgs system (as shown in Eq. ftL&)) 
satisfies a priori all the qualitative requirements (l)-(4), including the condition oJVckm 7^ 1- 1 
now discuss that this system works well even quantitatively. 

With the six effective Yukawa couplings shown in Eq. ([15]) , the Dirac mass matrices of quarks 
and leptons of the three families at the unification scale take the form: 



U 



( 



-e' 





mu, 



D 



I e' + 77' 

— e' + 7/ e + 77 
V -e + 77 1 



m D , 



N 



/ -3e' 

3e' 
\ 3e + a 




mu, 



/ - 3e' 
3e' + 77' 
V 3e + rj 



7/ 

-3e + r] I m D . 



(16) 



Here the matrices are multiplied by left-handed fermion fields from the left and by anti-fermion 
fields from the right. (U, D) stand for the mass matrices of up and down quarks, while (N, L) are 
the Dirac mass matrices of the neutrinos and the charged leptons. The entries 1, e,and a arise 
respectively from the ^33,023 and h 2 s terms in Eq. (|T5|), while 77 entering into D and L receives 
contributions from both g 2 3 and /i-23! thus f] 7^ a. Similarly 7/ and e' arise from g 12 and terms 
respectively. Note the quark-lepton correlations between U and iV as well as D and L arise because 
of SU(4) C , while the up-down correlations between U and D as well as iV and L arise because 
of SU(2)£xSU(2)#. Thus, these correlations emerge just because of the symmetry property of 
G(224). The relative factor of —3 between quarks and leptons involving the e entry reflects the 
fact that (45h) is proportional to (B-L), while the antisymmetry in this entry arises from the 
group structure of SO(10), as explained above 7 . As we will see, this e-entry helps to account for 
(a) the differences between m s and m M , (b) that between and m e , and most important, (c) 
the suppression of V c b together with the enhancement of the v^-Vt oscillation angle. 

The mass matrices in Eq. (|16D contain 7 parameters fj e, a, r/, m D = h 33 (10,^), mu = h 33 (10u), 



most likely to arise from an underlying 45 of SO(10) (rather than 210); in this case, the couplings of xh must be 
antisymmetric like that of 10 h ■ 45//. 

8 Of these, mj, w can in fact be estimated to within 20% accuracy by either using the argument of radiative 
electroweak symmetry breaking, or some promising string solutions (see e.g. Rcf. 
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rf and e'. These may be determined by using, for example, the following input values: m]? hys = 174 
GeV, m c (m c ) = 1.37 GeV, m s (l GeV) = 110-116 MeV @, m u (l GeV) « 6 MeV and the observed 
masses of e, \i and r, which lead to (see Ref. 0], for details): 



a 



0.110 , r? ~ 0.151 , e ~ - 0.095 , \r]'\ « 4.4 x 10" 3 and e' « 2 x 10' 
mc; ~ m t (Mu) ~ (100-120) GeV , m D ~ m 6 (M[/) ~ 1.5 GeV. 



(17) 

Here, I will assume, only for the sake of simplicity, as in Ref. |]14jl , that the parameters are 
real.0 Note that in accord with our general expectations discussed above, each of the parameters 
a, i] and e are found to be of order 1/10, as opposed to being 0(1) or O(10~ 2 ), compared to 
the leading (3,3)-element in Eq. (jlBf ). Having determined these parameters, we are led to a total 
of five predictions involving only the quarks (those for the leptons are listed separately): 



mi 



1 - 8e 2 ) ; thus m b (m b ) ~ (4.6-4.9) GeV 



w, 



cb\ 



\a — r)\ 



^m s /m b 



rj + e 



1/2 



- Jm c /m t 



rj — e 
m d (lGeV) ~ 8 MeV 



o + e 



a — e 



1/2 



0.045 



(19) 
(20) 



\]m d /m s - e lrt> ^m u /m c 



\V ub /V cb \ ~ \]m u /m c 



0.07. 



(21) 



[22) 



In making these predictions, we have extrapolated the GUT-scale values down to low energies 
using azimz) = 0.118, a SUSY threshold of 500 GeV and tan/3 = 5. The results depend weakly 
on these choices, assuming tan/5 w 2-30. Further, the Dirac masses and mixings of the neutrinos 
and the mixings of the charged leptons also get determined. We obtain: 



m°(Mu) w 100-120 GeV; (M v ) ~ 8 GeV 



D 



(23) 



9,.„ ~ — 3e + w Jm^/m T 



■fj,T 



3e + 77 



3e + 77 



1/2 



0.437 



[9e' 2 /(9e 2 - (t 2 )]!^^ ~ 0.4 MeV 



(24) 
(25) 



3e' 



77' + 3e' 



1/2 



nT'e/n^n — 0.85 Jm e /m 



0.06 



(26) 



9 Babu and I have recently studied supersymmetric CP violation within the G(224)/SO(10) framework, by 
using precisely the fermion mass- matrices as in Eq. ([ji]). We have observed |33) that complexification of the 
parameters can lead to observed CP violation, without upsetting in the least the success of Ref. ji3 (i.e. of the 
fermion mass-matrices of Eq. (|l6|)) in describing the masses and mixings of all fermions, including neutrinos. Even 
with complexification the relative signs and the approximate magnitudes of the real parts of the parameters must 
be the same as in Eq. (17), to retain the success. 

10 This is one characteristic difference between our work and that of Ref. ]43| , where the (2,3)-element is even 
bigger than the (3,3). 
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er 



0.85 



In evaluating 9 e efl , we have assumed e' and rf to be relatively positive. 

Given the bizarre pattern of quark and lepton masses and mixings, it seems remarkable that 
the simple and economical pattern of fermion mass-matrices, motivated in part by the assumption 
of flavor symmetries 6 which distinguish between the three families and in large part by the group 
theory of G(224)/SO(10), gives an overall fit to all of them [Eqs. fll~8| ) through fl2~2"|) l which is good 
to within 10%. This includes the two successful predictions on m b and V c b [Eqs.(|18"D and (|19|)1. 
Note that in supersymmetric unified theories, the "observed" value of mb(mb) and renormalization- 
group studies suggest that, for a wide range of the parameter tan/3, m' should in fact be about 
10-20% lower than mj ||45|| . This is neatly explained by the relation: « m°(l — 8e 2 ) [Eq. 
([18|)1, where exact equality holds in the limit e — > (due to SU(4)-color), while the decrease of 
m" compared to mj by 8e 2 ~ 10% is precisely because the off-diagonal e-entry is proportional to 



Specially intriguing is the result on V c b ~ 0.045 which compares well with the observed value 
of ~ 0.04. The suppression of V c b, compared to the value of 0.17 ± 0.06 obtained from Eq. flT3|), 
is now possible because the mass matrices [Eq. (|16D] contain an antisymmetric component oc e. 



That corrects the square-root formula 6 s b = ym s /m b [appropriate for symmetric matrices, see Eq. 

flTID I by the asymmetry factor | (77-he) / (77 — e)) 1 / 2 [see Eq. (19)], and similarly for the angle 9 ct . This 
factor suppresses V c b if r\ and e have opposite signs. The interesting point is that, the same feature 
necessarily enhances the corresponding mixing angle Q in the leptonic sector, since the asymmetry 
factor in this case is given by [(— 3e + 77) / (3e + i])] 1 ^ 2 [see Eq. (24)]. This enhancement of 8 e 
helps to account for the nearly maximal oscillation angle observed at SuperK (as discussed below). 
This intriguing correlation between the mixing angles in the quark versus leptonic sectors — that is 
suppression of one implying enhancement of the other — has become possible only because of the 
e-contribution, which is simultaneously antisymmetric and is proportional to B L. That in turn 
becomes possible because of the group-property of SO(10) or a string-derived G(224) 7 . 

Taking stock, we see an impressive set of facts in favor of having B-L as a gauge symmetry 
and in fact for the full SU(4)-color-symmetry. These include: (i) the suppression of V c b, together 
with the enhancement of 9 , mentioned above; (ii) the successful relation ~ mP T (l — 8e 2 ); 
(iii) the usefulness again of the SU(4)-color-relation ^(^Dirac) ~ m t i n accounting for m(z/£) (see 
Section 4); (iv) the agreement of the relation \rn? s /rn?\ = |(e 2 — r] 2 )/(9e 2 — rf)\ with the data, in 
that the ratio is naturally less than 1, if n ~ e [The presence of 9e 2 in the denominator is because 
the off-diagonal entry is proportional to B-L.]; and finally (v), the need for (B-L) — as a local 
symmetry, to implement baryogenesis via leptogenesis, as noted in Section 1. 

Turning to neutrino masses, while all the entries in the Dirac mass matrix iV are now fixed, 
to obtain the parameters for the light neutrinos, one needs to specify those of the Majorana mass 
matrix of the RH neutrinos (^ ,M ' T ). Guided by economy and the assumption of hierarchy, we 
consider the following pattern ||14|| : 



As discussed in Section |, the magnitude of M R rs (5-10) x 10 14 GeV can quite plausibly be 
justified in the context of supersymmetric unification 5 [e.g. by using M « M st m 4 x 10 17 GeV in 



B-L [see Eq. (0)]. 





/ x z 
y 
\ z y 1 




(28) 
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Eq. (H)]. To the same extent, the magnitude of m{y T ) ~ (1/10-1/30) eV, which is consistent with 



the SuperK value, can also be anticipated by allowing for — v T mixing [see Ref. [14j]]. Thus 
there are effectively three new parameters: x, y, and z. Since there are six observables for the 
three light neutrinos, one can expect three predictions. These may be taken to be 9^ c Ut , rn VT [see 
Eq. (|i~0D1, and for example 9° s e l ■ 

Assuming successively hierarchical entries as for the Dirac mass matrices, we presume that 
\y\ ~ 1/10, \z\ < \y\/10 and \x\ < z 2 . Now given that m(v T ) ~ 1/20 eV [as estimated in Eq. 



([T0|)1, the MSW solution for the solar neutrino puzzle jKJ suggests that m{v^)/m{i/ T ) ~ 1/8-1/20. 
With hierarchical neutrino masses, the higher value of the mass-ratio (like 1/8) holds only for the 
large angle MSW solution (see below). With the mass-ratio being in the range of 1/8-1/20, one 
obtains: \y\ « (1/17 to 1/21), with y having the same sign as e [see Eq. (0)]. This solution for 
y obtains only by assuming that y has a hierarchical value 0(1/10) rather than 0(1). Combining 
now with the mixing in the fi-r sector determined above [see Eq. (|4])], one can then determine 
the Vp-Vr oscillation angle. The two predictions of the model for the neutrino-system are then: 

m(v T ) « (1/10- 1/30) eV (29) 

9™ ^9^-9^^ (o-437 + ^) • (30) 

Thus, 

sin 2 20°^ Vt = (0.99, 0.975, 0.92, 0.87) (31) 

for 

mJm U3 = (1/8, 1/10, 1/15, 1/20) . (32) 

Both of these predictions are extremely successful .0 

Note the interesting point that the MSW solution, and the requirement that \y\ should have a 
natural hierarchical value (as mentioned above), lead to y having the same sign as e. Now, that (it 
turns out) implies that the two contributions in Eq. (|30|) must add rather than subtract, leading 



to an almost maximal oscillation angle |TJj]. The other factor contributing to the enhancement 
of 9° bc Vt is, of course, also the asymmetry-ratio which increases |6^ T | from 0.25 to 0.437 [see Eq. 
(124D1- We see that one can derive rather plausibly a large v^-v T oscillation angle sin 2 29° sc Ut > 0.92, 



together with an understanding of hierarchical masses and mixings of the quarks and the charged 
leptons, while maintaining a large hierarchy in the seesaw derived neutrino masses (m U2 /m U3 = 
1/8-1/15), all within a unified framework including both quarks and leptons. In the example 
exhibited here, the mixing angles for the mass eigenstates of neither the neutrinos nor the charged 
leptons are really large, in that 0j T ~ 0.437 ~ 23° and 0£ r ~ (0.22-0.35) w (13-20.5)°, yet the 
oscillation angle obtained by combining the two is near-maximal. This contrasts with most works 
in the literature in which a large oscillation angle is obtained either entirely from the neutrino 
sector (with nearly degenerate neutrinos) or almost entirely from the charged lepton sector. 

Small Versus Large Angle MSW Solutions 

In considerations of z/ e -z/^ and v e -v T oscillation angles, tiny intrinsic non-diagonal Majorana masses 
~ 10~ 3 eV of the LH neutrinos leading to v^v^ and i/£i/£-mixings, which can far exceed those 

11 In writing Eq. (JIl|), the small angle approximation exhibited in Eq. ( ^0| ) is replaced by the more precise 
expression, given in Eq. (12) of Ref. |14|, with the further understanding that ^Jm^/m T appearing in Eq. (12) 



(of Ref. UM) is replaced by the /i-r mixing angle w 0.437. 
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induced by the standard see-saw mechanism, can be rather important, especially for z/ e -z/ M mixing. 
As explained below, such intrinsic masses can arise quite naturally through higher dimensional 
operators and can lead to the large angle MSW solution of the solar neutrino puzzle. 

Let us first ignore the intrinsic Majorana masses of the LH neutrinos and include only those 
that arise through the standard see-saw mechanism, involving the superheavy Majorana masses 
of the RH neutrinos, with a pattern given, for example, by Eq. (f28|). Note that, while Mr « 
(5-15) x 10 14 GeV and y « —1/20 are better determined, the parameters x and z can not be 
obtained reliably at present because very little is known about observables involving v e . Taking, 
for concreteness, m Ue « (10~ 5 -10~ 4 ) (1 to few)) eV and 9°* c « 6 e eT — 6 u eT « 10~ 3 ± 0.03 as inputs, 
we obtain: z ~ (1-5) x 10~ 3 and x ~ (1 to few)(10~ 6 -10~ 5 ), in accord with the guidelines of 
|*| ~ \y\/10 and \x\ ~ z 2 . This in turn yields: 9™ « 0^ - 9^ « 0.06 ± 0.015. Note that the 
mass of m v ~ 3 x 10~ 3 eV, that follows from a natural hierarchical value for y ~ —(1/20), and 
8 efJi as above, go well with the small angle MSW explanation of the solar neutrino puzzle. In 
short the framework presented so far, that neglects intrinsic Majorana masses of the LH neutrinos 
altogether, generically tends to yield the small angle MSW solution. 

As alluded to above, we now observe that small intrinsic non-seesaw masses of the LH neutrinos 
~ 10 -3 eV, which could mix u e L and v^l, can, however, arise quite naturally through higher di- 
mensional operators in the superpotential of the form[^ W D fti2l6il62l6# 16#10#10h/M gut . 
One can verify that such a term would lead to an intrinsic Majorana mixing mass term of 
the form m^u^i/^, with a strength given by m^} m k 12 ((16#)/M GU t) 2 (175 GeV) 2 /M GU T ~ 
(1.5-6) x 10- 3 eV, where we have put (16 H ) « (1-2)M GU t and M GU t ~ 2 x 10 16 GeV. Such an 
intrinsic Majorana mixing mass ~ 10~ 3 eV, though small, is still much larger than what one would 
get for the corresponding term from the standard see-saw mechanism. Now, as discussed above, 
the diagonal mass-term, arising from the standard see-saw mechanism can naturally be of 

order (3-8)xl0 -3 eV (for \y\ m 1/20 to 1/15, say). In addition, the intrinsic contribution of the 
type mentioned above may in general also contribute to the diagonal {vl, v l) mass (depending upon 
flavor symmetries) which can be (few)xlO -3 eV. Thus, taking the net values of 1TI22 ~ (6-7) x 10~ 3 
eV (say), m$ « (3-4) x 10" 3 eV, and mf^ < (1-2) x 10 3 eV, which are all very plausible, we 
obtain m v ~ (6-7) x 10 -3 eV, m Ve ~ 1 x 10~ 3 eV, so that Amf 2 « (3.6-5) x 10~ 5 eV 2 , and 
sin 2 26*°2 C ~ 0.6-0.7. This goes well with the large angle MSW solution of the solar neutrino 
puzzle, which is now favored over the small angle solution by the SuperK data fT? . 



In summary, the intrinsic non-seesaw contribution to the Majorana masses of the LH neutrinos 
quite plausibly has the right magnitude for v^v^ mixing, so as to lead to the rather large oscillation 
angle as mentioned above, in accord with the data. In contrast to the case of the v^-v^ oscillation 
angle, however, given the smallness of the entries involving the first two families, the relatively large 
angle solution for v e —v^ oscillation may not be regarded as a firm prediction of the SO(10)/G(224)- 
framework presented here. It is nevertheless a very reasonable possibility. 

It is worth noting that although the superheavy Majorana masses of the RH neutrinos can- 
not be observed directly, they can be of cosmological significance. The pattern given above 
and the arguments given in Section ||] and in this section suggests that M{y T ^) w (5-15) x 10 14 
GeV, M(z/£) w (1-4) x 10 12 GeV (for \y\ w 1/20); and M(u e R ) ~ (1/2-10) x 10 9 GeV (for 
x ~ (l/2-10)10~ 6 > z 2 ). A mass of v e R ~ 10 9 GeV is of the right magnitude for producing v e R 
following reheating and inducing lepton asymmetry in v e R decay into H° + 1/^, that is subsequently 



12 Such a term can be induced in the presence of, for example, a singlet S and a ten-plet (denoted by 10), both 
having GUT-scale masses, and possessing renormalizable couplings of the form ail6il6//10, blOlOnS, MsSS 

^GUT 



and M10 2 . In this case, n 12 /M^ VT = a 1 a 2 b/(M 2 M s ). 
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converted into baryon asymmetry by the electroweak sphalerons |]T6|, [L7| . 

In summary, we have proposed an economical and predictive pattern for the Dirac mass ma- 
trices, within the SO(10)/G(224)-framework, which is remarkably successful in describing the 
observed masses and mixings of all the quarks and charged leptons. It leads to five predictions 
for just the quark- system, all of which agree with observation to within 10%. The same pat- 
tern, supplemented with a similar structure for the Majorana mass matrix, accounts for both the 
nearly-maximal v^-v T oscillation angle and a (mass) 2 -difference Am 2 (^i/ T ) ~ (1/20 eV) 2 , sug- 
gested by the SuperK data. Given this degree of success, it makes good sense to study proton 
decay concretely within this SO(10)/G(224)-framework. The results of this study |TJ|, |TB|] are 
presented in the next section, together with an update. 

Before turning to proton decay, it is worth noting that much of our discussion of fermion 
masses and mixings, including those of the neutrinos, is essentially unaltered if we go to the limit 
e' — > of Eq. (28). This limit clearly involves: 



0, 9 C ~ ^m d /m s , m Ve = 0, 9% = 6 u eT = 



\ Vut> \ ~ y \/ m d/m b (m 8 /m b ) ~ (2.1)(0.039)(0.023) ~ 0.0019 . (33) 

All other predictions remain unaltered. Now, among the observed quantities in the list above, 
Q c — \Jrn d /m s is a good result. Considering that m u /m t « 10~ 5 , m u = is also a pretty 
good result. There are of course plausible small corrections which could arise through Planck 
scale physics; these could induce a small value for m u through the (l,l)-entry 5 ~ 10 -5 . For 
considerations of proton decay, it is worth distinguishing between these two extreme variants 
which we will refer to I and II respectively. 

Case I : e' « 2 x 10" 4 ,6 = 

Case II: 5 « 10~ 5 , e' = . (34) 

It is worth noting that the observed value of \V u b\ ~ 0.003 favors a non-zero value of e' (~ (1- 
2) x 10~ 4 ). Thus, in reality, e' may not be zero, but it may lie in between the two extreme values 
listed above. In this case, the predicted proton lifetime for the standard d = 5 operators would 
be intermediate between those for the two cases, presented in Section ||. 

6 Expectations for Proton Decay in Supersymmetric Unified 
Theories 

6.1 Preliminaries 

Turning to the main purpose of this talk, I present now the reason why the unification framework 
based on SUSY SO(10) or G(224), together with the understanding of fermion masses and mixings 
discussed above, strongly suggest that proton decay should be imminent. 

Recall that supersymmetric unified theories (GUTs) introduce two new features to proton 
decay: (i) First, by raising Mx to a higher value of about 2 x 10 16 GeV (contrast with the non- 
supersymmetric case of nearly 3 x 10 14 GeV), they strongly suppress the gauge-boson-mediated 
d = 6 proton decay operators, for which e + 7r° would have been the dominant mode (for this 
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case, one typically obtains: T l (p — > e + 7r°)|d =6 « 10 35±1 years), (ii) Second, they generate d = 5 



proton decay operators ]T9) of the form QiQjQkQi/M in the superpotential, through the exchange 
of color triplet Higginos, which are the GUT partners of the standard Higgs(ino) doublets, such as 
those in the 5 + 5 of SU(5) or the 10 of SO(10). Assuming that a suitable doublet-triplet splitting 
mechanism provides heavy GUT-scale masses to these color triplets and at the same time light 
masses to the doublets (see e.g, the Appendix), these "standard" d — 5 operators, suppressed 
by just one power of the heavy mass and the small Yukawa couplings, are found to provide the 



dominant mechanism for proton decay in supersymmetric GUT |48|. [49|. |50], |5T], p2 | . 

Now, owing to (a) Bose symmetry of the superfields in QQQL/M, (b) color antisymmetry, 
and especially (c) the hierarchical Yukawa couplings of the Higgs doublets, it turns out that these 
standard d = 5 operators lead to dominant VK + and comparable Vn + modes, but in all cases to 
highly suppressed e + 7r°, e + K° and even fi + K° modes. For instance, for minimal SUSY SU(5), 
one obtains (with tan/3 < 20, say): 

[T{^K°)/T{vK + )] s s ^ ~ K/(m c sin 2 9 C )] 2 R » KT 3 , (35) 

where R ~ 0.1 is the ratio of the relevant |matrix element| 2 x (phase space), for the two modes. 

It was recently pointed out that in SUSY unified theories based on SO(10) or G(224), which 
assign heavy Majorana masses to the RH neutrinos, there exists a new set of color triplets and 



thereby very likely a new source of d = 5 proton decay operators [2C]. For instance, in the con- 
text of the minimal set of Higgs multiplets {45//, 16//, 16// and 10//} (see Section |5]), these 
new d = 5 operators arise by combining three effective couplings introduced before: — i.e., (a) the 
couplings /jjl6jl6jl6#16#/M [see Eq. ([?])] that are required to assign Majorana masses to the 
RH neutrinos, (b) the couplings gj J 16jl6 J 16//16///M, which are needed to generate non-trivial 
CKM mixings [see Eq. ( fT5|) 1, and (c) the mass term M 16 16#16#. For the fij couplings, there are 
two possible SO(10)-contractions (leading to a 45 or a 1) for the pair 16,16//, both of which con- 
tribute to the Majorana masses of the RH neutrinos, but only the non-singlet contraction (leading 
to 45), would contribute to d — 5 proton decay operator. In the presence of non-perturbative 
quantum gravity, one would in general expect the two contractions to have comparable strength. 
Furthermore, the couplings of 45's lying in the string-tower or possibly below the string-scale, and 
likewise of singlets, to the 16, ■ 16//-pair, would respectively generate the two contractions. It thus 
seems most likely that both contractions would be present, having comparable strength. Allowing 
for a difference between the relevant projection factors for vr masses versus proton decay, and 
also for the fact that both contractions contribute to the former, but only the non-singlet one (i.e. 
45) to the latter, we would set the relevant fij coupling for proton decay to be (fij) p = (fij)v • K, 
where (fij) u defined in Section f| directly yields Ur - masses [see Eq. (§)]; and K is a relative 
factor, which generically is expected to be of order unity. As a plausible range, we will take 
X « 1/5 to 2 (say). In the presence of the non-singlet contraction, the color-triplet Higginos in 
16h and 16// of mass M 16 can be exchanged between and q^-pairs (correspondingly, for 
G(224), the color triplets would arise from (1,2,4)// and (1,2,4)//). This exchange generates a 
new set of d = 5 operators in the superpotential of the form 

W new oc (f tJ ) u g kl K (16, 16,) (16, 16,) (16//) (16 H }/M 2 x (1/M 16 ), (36) 



13 The origin of the new d = 5 operators in the context of other Higgs multiplets, in particular in the cases where 
126// and 126// are used to break B-L, has been discussed in Ref. p0[ . 

1 For the special case of K = (which would arise if only the singlet-contraction of (16j ■ 16//) contributes), the 
new d — 5 operators shown in Eq. (pq) would not, of course, contribute to proton decay. 
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which induce proton decay. Note that these operators depend, through the couplings and g^i, 
both on the Majorana and on the Dirac masses of the respective fermions. This is why within 
SUSY SO(10) or G(224), if the generic case of K ^ holds, proton decay gets intimately linked 
to the masses and mixings of all fermions, including neutrinos. 



6.2 Framework for Calculating Proton Decay Rate 



To establish notations, consider the case of minimal SUSY SU(5) and, as an example, the process 
cd — > s&p, which induces p — > V fl K + . Let the strength of the corresponding d = 5 operator, 
multiplied by the product of the CKM mixing elements entering into wino-exchange vertices, 
(which in this case is sin Q c cos 9c) be denoted by A. Thus (putting cos 9c = 1), one obtains: 



A Bd (SU(5)) = (h u 22 h d 12 /M Hc ) sm6 c 



(m c m s sm 2 9 c /vl) (t&n(3/M Hc ) 
(1.9 x 1(T 8 ) (t&nf3/M Hc ) 
(2 x 10~ 24 GeV- 1 ) (tan/3/2) (2 x 10 16 GeV/M, 



(37) 



where tan/? = v u /vd, and we have put v u = 174 GeV and the fermion masses extrapolated to 
the unification-scale — i.e. m c ~ 300 MeV and m s ~ 40 MeV. The amplitude for the associated 
four-fermion process dus — > 



is given by: 
A 5 (dus -> 77 M ) = x (2/) 



(3? 



where / is the loop-factor associated with wino-dressing. Assuming <C rriq ~ mj, one gets: 
/ ~ (m^/m?) (a 2 /47r). Using the amplitude for (du)(sui), as in Eq. fl38D , (£ = p or r), and 
the recently obtained matrix element and renormalization effects (see below), one then obtains 
m 151], H 0, 0: 



r _1 (p -> z7 T K + ) sa (0.15 x 10 31 ) years x (0.32/ A L ) 



x 



0.93 
^47 



0.014 GeV 



31 2 



(1/6) 



771; 



Ll.2TeV 



(39) 



2 x 10~ 24 GeV 



7—1 



A(I7) 



Here (3h denotes the hadronic matrix element defined by (3nUL(k) = e a /j 7 ^0|(d2^)wX|p, fc^. While 
the range (3jj = (0.003-0.03) GeV 3 has been used in the past ]5T| , given that one lattice calculation 
yields (3h = (5.6 ± 0.5) x 10 -3 GeV 3 |kJ, and a recent improved calculation yields (3h ~ 0.014 
GeV 3 |54| (whose systematic errors that may arise from scaling violations and quenching are hard 
to estimate [0), we will take as a conservative, but plausible, range for (3h to be given by (0.014 
GeV 3 )(l/2 - 2). (Compare this with the range for f3 H = (0.006 GeV 3 )(l/2 - 2) as used in Ref. 
P^j). As denotes the short- distance renormalization effect for the d — 5 operator which arises 
owing to extrapolation between the GUT and the SUSY-breaking scales [[49], 51], [55| . The average 
value of As = 0.67, given in Ref. [5"I| for m t = 100 GeV, has been used in most early estimates. 
For m t = 175 GeV, one would, however, have As ~ 0.93 to 1.2 []55| . Conservatively, I would use 
As = 0.93; this would enhance the rate by a factor of two compared with previous estimates. Al 
denotes the long-distance renormalization effect of the d = 6 operator due to QCD interaction that 



arises due to extrapolation between the SUSY breaking scale and 1 GeV j4~9| . Using the two- loop 
expression for Ai together with the two-loop value for a 3 , Babu and I find: Al ~ 0.32, in 



21 



contrast to Al « 0.22, used in previous worksQ. In what follows, I would use Al m 0.32. This 
by itself would also increase the rate by a factor of (0.32/0.22) 2 «2, compared to the previous 
estimates fl49|, pOl pTL [521 [14], p75| . Including the enhancements in both As and Al, we thus see 



that the net increase in the proton decay rate solely due to new evaluation of renormalization 
effects is nearly a factor of four, compared to the previous estimates (including that in Ref. ||14|| ). 

Note that the familiar factors that appear in the expression for proton lifetime — i.e., Mh c , 
(1 + y tc ) representing the interference between the t and c contributions, and tan/? (see e.g. Ref. 
plf and discussion in the Appendix of Ref. |14|) — are all effectively contained in A(v). In Ref. 



14j| , guided by the demand of naturalness (i.e. absence of excessive fine tuning) in obtaining the 
Higgs boson mass, squark masses were assumed to lie in the range of 1 TeV(l/ v / 2 - V2), so that 
rriq < 1.4TeV. Recent work, based on the notion of focus point supersymmetry however suggests 
that squarks may be considerably heavier without conflicting with the demands of naturalness 

58| . In the interest of obtaining a conservative upper limit on proton lifetime, we will therefore 



allow squark masses to be as heavy as about 2.4 TeV and as light as perhaps 600 GeV. 

Allowing for plausible and rather generous uncertainties in the matrix element and the spec- 
trum we take: 

fa = (0.014 GeV 3 ) (1/2-2) 
ma/niq = 1/6(1/2-2), and m~ q w m i « 1.2 TeV (1/2 - 2) . (40) 
Using Eqs. (39-40), we get: 

F~ l (p -> V T K + ) w (0.15 x 10 31 years) [2 x 1(T 24 GeV _1 /^fa) ] 2 x {64-1/64}. (41) 



Note that the curly bracket would acquire its upper-end value of 64, which would serve towards 
maximizing proton lifetime, only provided all the uncertainties in Eq. ([41]) are stretched to the 
extreme so that flu = 0.007 GeV 3 , m^/rriq m 1/12 and rriq ps 2.4 TeV. This relation, as well as 
Eq. (39) are general, depending only on A(vi) and on the range of parameters given in Eq. (40). 
They can thus be used for both SU(5) and SO(10). 

15 In most previous works starting with Ref. |49) through p^j, as well as in Refs. |Q and |l8| , the one-loop 
value of Al was taken to be 0.22. It was, however, noted in Refs. |5(| and |57j that there is a numerical error in 
the evaluation of the one-loop expression for Al p^ |, and that the correct value for ^(one — loop) 0.43 (this 
remained unnoticed by most authors). The two-loop value for Al (as stated above) is nearly 0.32, which is lower 
than 0.43 but higher that the previously used value of 0.22. 

16 We remark that if the recently reported (g-2)-anomaly for the muon |5{| , together with reevaluation of the 
contribution from light by light-scattering IHj], is attributed to supersymmetry [Ell , one would need to have 
extremely light s-fermions [i.e. mj ss 200 - 400 GeV (say) and correspondingly, for promising mechanisms of 
SUSY- breaking, rriq < 300 — 600 GeV (say)], and simultaneously relatively large tan/3(ss 6-24). However, not 
worrying about grand unification, such light s-fermions, together with large or very large tan/3 would typically be 
in gross conflict with the limits on the edm's of the neutron and the electron, unless one can explain naturally 
the occurrence of minuscule phases 1/200 to 1/500) and/or large cancellation. Thus, if the (g — 2) At -anomaly 
turns out to be real, it may well find a non-supersymmetric explanation, in accord with the edm-constraints which 
ordinarily seem to suggest that squarks are (at least) moderately heavy (rriq ^ 0.6 — 1 TeV, say), and tan/3 is 
not too large (< 3 to 10, say). We mention in passing that the extra vector — like matter — specially a 16 + 16 of 
SO(10) — as proposed in the so-called extended supersymmetric standard model (ESSM) |H], |6^], with the heavy 
lepton mass being of order 200 GeV, can provide such an explanation Q. Motivations for the case of ESSM, 
based on the need for (a) removing the mismatch between MSSM and string unification scales, and (b) dilaton- 
stabilization, have been noted in Ref. pi] . Since ESSM is an interesting and viable variant of MSSM, and would 
have important implications for proton decay, we will present the results for expected proton decay rates for the 
cases of both MSSM and ESSM in the discussion to follow. 
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The experimental lower limit on the inverse rate for the vK + modes is given by Ref. 164 



£ T(p - V t K- 



> 1.9 x 10 33 years. 



(42) 



expt 



Allowing for all the uncertainties to stretch in the same direction (in this case, the curly bracket 
= 64), and assuming that just one neutrino flavor (e.g. for SU(5)) dominates, the observed 
limit (Eq. (f|2|)) provides an upper bound on the amplitude^]: 



A(v e ) < 0.46 x KT 24 GeV -1 



(43) 



which holds for both SU(5) and SO(10). Recent theoretical analyses based on LEP-limit on Higgs 
mass (> 114 GeV), together with certain assumptions about MSSM parameters (as in CMSSM) 
and/or constraint from muon g-2 anomaly [59] suggest that tan/? > 3 to 5 |SS |. In the interest 
of getting a conservative upper limit on proton lifetime, we will therefore conservative 
lower limit, tan/3 > 3. We will however exhibit relevant results often as a function of tan/? and 
exhibit proton lifetimes corresponding to higher values of tan/? as well. For minimal SU(5), using 
Eqs. (37) and (|43"D and, conservatively tan/3 > 3, one obtains a lower limit on M HC given by: 



M HC > 13 x 10 ib GeV (SUSY SU(5)) 



(44) 



At the same time, gauge coupling unification in SUSY SU(5) strongly suggests Mhc < (1/2-1) x 
10 16 GeV. (See Ref. [66] where an even more stringent upper bound on Mhc is suggested.) Thus 
we already see a conflict, in the case of minimal SUSY SU(5), between the experimental limit on 
proton lifetime on the one hand, and coupling unification and constraint on tan/? on the other 
hand. To see this conflict another way, if we keep Mhc < 10 16 GeV (for the sake of coupling 
unification) we obtain from Eq. (37): A(SU(5)) > 5.7 x KT 24 GeV -1 (tan/?/3). Using Eq. ©, 
this in turn implies that 



r -1 (p -> VK + ) < 1.2 x 10 31 years x (3/ tan/?) 2 (SUSY SU(5)) 



(45) 



For tan/? > 3, a lifetime of 1.2 x 10 31 years is thus a most conservative upper limit. In practice, 
it is unlikely that all the uncertainties, including these in Mhc an d tan/?, would stretch in the 
same direction to nearly extreme values so as to prolong proton lifetime. Given the experimental 
lower limit [Eq. fl4"2|)l, we see that minimal SUSY SU(5) is already excluded by a large margin by 
proton decay-searches. This is in full accord with the conclusion reached by other authors (see 
especially Ref. [p6]]). We have of course noted in Section [| that SUSY SU(5) does not go well 
with neutrino oscillations observed at SuperK. 

Now, to discuss proton decay in the context of supersymmetric SO (10), it is necessary to 
discuss first the mechanism for doublet-triplet splitting. Details of this discussion may be found 
in Ref. 



TM. A synopsis is presented in the Appendix. 



6.3 Proton Decay in Supersymmetric SO(10) 

The calculation of the amplitudes A st< i and A new for the standard and the new operators for 
the SO(10) model, are given in detail in Ref. [14]. Here, I will present only the results. It is 
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If there are sub-dominant ViK + modes with branching ratio R, the right side of Eq. (E3I) should be divided 



by y/T+R. 
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found that the four amplitudes A st d(v T K + ), A st d(T> fl K + ), A ncw (v T K + ) and A new (v^K + ) are in 
fact very comparable to each other, within about a factor of two to five, either way. Since there 
is no reason to expect a near cancellation between the standard and the new operators, especially 
for both V T K + and V fl K + modes, we expect the net amplitude (standard + new) to be in the 



range exhibited by either one. Following Ref. [I4]j, I therefore present the contributions from the 
standard and the new operators separately. 

One important consequence of the doublet-triplet splitting mechanism for SO(10) outlined 
briefly in the appendix and in more detail in Ref. |L4]| is that the standard d = 5 proton decay 



operators become inversely proportional to M c s = [A(45#)] 2 / M W r ~ M x /M\qi, rather than 
to Mh c - Here, M\qi represents the mass of 10'^, that enters into the D-T splitting mechanism 
through effective coupling \1Q H 45 H 10' H in the superpotential [see Appendix, Eq. ([Alp ]. As noted 
in Ref. [fL4}| , M 1( y can be naturally suppressed (due to flavor symmetries) compared to M x , and 
thus M e g correspondingly larger than Mx by even one to three orders of magnitude. It should 
be stressed that M e fj does not represent the physical masses of the color triplets or of the other 
particles in the theory. It is simply a parameter of order M x /Miq/. Thus values of M e R, close 
to or even exceeding the Planck scale, do not in any way imply large corrections from quantum 
gravity. Now accompanying the suppression due to M eff , the standard proton decay amplitudes 
for SO(10) possess an intrinsic enhancement as well, compared to those for SU(5), owing primarily 
due to differences in their Yukawa couplings for the up sector (see Appendix C of Ref. fll4jl ). As a 



result of this enhancement, combined with the suppression due to higher values of M e g, a typical 
standard d = 5 amplitude for SO(10) is given by (see Appendix C of Ref. flip 

^(^ + )S (10) ^(^33/M cff )(2xlO- 5 ), 

which should be compared with A^K+ff/^ « (1.9 x l(T 8 )(tan/?/^fec) [ see E Q- (37)]. Note, 
taking h^ 3 ~ 1/4, the ratio of a typical SO(10) over SU(5) amplitude is given by 
(M# c /M eff )(88)(3/ tan/3). Thus the enhancement by a factor of about 88 (for tan/3 = 3), of 
the SO (10) compared to the SU(5) amplitude, is compensated in part by the suppression that 
arises from M e g being larger than M# c . 

In addition, note that in contrast to the case of SU(5), the SO(10) amplitude does not depend 
explicitly on tan/3. The reason is this: if the fermions acquire masses only through the 10 h in 
SO (10), as is well known, the up and down quark Yukawa couplings will be equal. By itself, 
it would lead to a large value of tan/3 = m t /mb ~ 60 and thereby to a large enhancement in 
proton decay amplitude. Furthermore, it would also lead to the bad relations: m c /m s = m t /mb 
and Vckm = 1- However, in the presence of additional Higgs multiplets, in particular with the 
mixing of (16^-)^ with 10# (see Appendix and Section (a) tan/3 can get lowered to values like 
3-20, (b) fermion masses get contributions from both (16#)d and (10#), which correct all the bad 
relations stated above, and simultaneously (c) the explicit dependence of A on tan/3 disappears. 
It reappears, however, through restriction on threshold corrections, discussed below. 

Although M e g can far exceed Mx, it still gets bounded from above by demanding that coupling 
unification, as observed PI, should emerge as a natural prediction of the theory as opposed to 



18 For instance, in the absence of GUT-scale threshold corrections, the MSSM value of az(rnz)MSSM , assuming 
coupling unification, is given by ct^{rnz)° MSSM — 0.125 ± 0.13 Q, which is about 5-8% higher than the observed 
value: a3(mz)° M sSM = 0.118±0.003 |0. We demand that this discrepancy should be accounted for accurately by 



a net negative contribution from D-T splitting and from "other" threshold corrections [see Appendix, Eq. (A4)], 
without involving large cancellations. That in fact does happen for the minimal Higgs system (45, 16, 16) (see Ref. 

ED- 
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being fortuitous. That in turn requires that there be no large (unpredicted) cancellation between 
GUT-scale threshold corrections to the gauge couplings that arise from splittings within different 
multiplets as well as from Planck scale physics. Following this point of view, we have argued (see 
Appendix) that the net "other" threshold corrections to Oi${mz) arising from the Higgs (in our 
case 45//, 16// and 16//) and the gauge multiplets should be negative, but conservatively and 
quite plausibly no more than about 10%, at the electroweak scale. This in turn restricts how big 
can be the threshold corrections to a^ivriz) that arise from (D-T) splitting (which is positive). 
Since the latter is proportional to ln(M e ff cosj/Mx) (see Appendix), we thus obtain an upper limit 
on M e g cos 7. For the simplest model of D-T splitting presented in Ref. |l|] and in the Appendix 
[Eq. (|A1|) 1, one obtains: cos 7 ~ (tan /3) / \m t / m b ) . An upper limit on M efT cos7 thus provides 



an upper limit on M e g which is inversely proportional to tan/3. In short, our demand of natural 
coupling unification, together with the simplest model of D-T splitting, introduces an implicit 
dependence on tan/3 into the lower limit of the SO(10)-amplitude — i.e. A(SO (10)) oc 1/M e g > 
[(a quantity) oc tan/3]. These considerations are reflected in the results given below. 

Assuming tan/3 > 3 and accurate coupling unification (as described above), one obtains for 
the case of MSSM, a conservative upper limit on M e g < 2.7 x 10 18 GeV (3/ tan/3) (see Appendix 
and Ref. |14]]). Using this upper limit, we obtain a lower limit for the standard proton decay 
amplitude given by 



A(V T K' 



Ifitd 



> 



(7.8 x KT^GeV- 1 
(3.3 x lO-^GeV -1 ; 



) (1/6-1/4) easel 
(1/6-1/2) case II 



SO(10)/MSSM, with 
tan/3 > 3 



(46) 



Substituting into Eq. fl4T| ) and adding the contribution from the second competing mode f^K* , 
with a typical branching ratio R ~ 0.3, we obtain 



Y-\vK + ) std < 



(0.18 x 10 31 years) (1.6-0.7) 
(0.4 x 10 31 years) (4-0.44) 



{64-1/64} 



SO(10)/MSSM, with 
tan/3 > 3 



(47) 



The upper and lower entries in Eqs. ( |46l) and ([|7]) correspond to the cases I and II of the fermion 
mass-matrix with the extreme values of e' — i.e. e' = 2x 10 -4 and e' = — respectively, (see Eq. 
(|34]) ). The uncertainty shown inside the square brackets correspond to that in the relative phases 



of the different contributions. The uncertainty of {64 to 1/64} arises from that in /3//, 



[m 



and rriq [see Eq. (0)]. Thus we find that for MSSM embedded in SO(10), for the two extreme 



values of e' 
satisfy: 



(cases I and II) as mentioned above, the inverse partial proton decay rate should 



Istd 



< 



< 



q-l +2.0 

0.20 x 10 31 -i-7 



0.32 x 10 



31 



+2.4 



years 



years 



0.2 x 10 33 years 
1 x 10 33 years 



SO(10)/MSSM, with 
tan/3 > 3 



(4f 



The central value of the upper limit in Eq. (|48| ) corresponds to taking the upper limit on M e g < 
2.7 x 10 18 GeV, which is obtained by restricting threshold corrections as described above (and in 
the Appendix) and by setting (conservatively) tan/3 > 3. The uncertainties of matrix element, 
spectrum and choice of phases are reflected in the exponents. The uncertainty in the most sensitive 
entry of the fermion mass matrix — i.e. e' — is incorporated (as regards obtaining an upper limit 
on the lifetime) by going from case I (with e' = 2 x 10~ 4 ) to case II (e' = 0). Note that this 
increases the lifetime by almost a factor of six. Any non- vanishing intermediate value of e' would 
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only shorten the lifetime compared to case II. In this sense, the larger of the two upper limits 
quoted above is rather conservative. We see that the predicted upper limit for case I of MSSM 
(with the extreme value of e' = 2 x 1CT 4 ) is lower than the empirical lower limit [Eq. fl43|)1 by a 
factor of ten, while that for case II, i.e. e' = (with all the uncertainties stretched as mentioned 
above) is about two times lower than the empirical lower limit. 

Thus the case of MSSM embedded in SO(10) is already tightly constrained, to the point of 
being disfavored, by the limit on proton lifetime. The constraint is of course augmented especially 
by our requirement of natural coupling unification which prohibits accidental large cancellation 
between different threshold corrections^ (see Appendix); and it will be even more severe, especially 
within the simplest mechanism of D-T splitting (as discussed in the Appendix), if tan (3 turns out 
to be larger than 5 (say). On the positive side, improvement in the current limit by a factor of 
even 2 to 3 ought to reveal proton decay, otherwise the case of MSSM embedded in SO(10), would 
be clearly excluded. 



6.4 The case of ESSM 

Before discussing the contribution of the new d — 5 operators to proton decay, an interesting 
possibility, mentioned in the introduction (and in footnote 16), that would be especially relevant 
in the context of proton decay, if tan/3 is large, is worth noting. This is the case of the extended 
supersymmetric standard model (ESSM), which introduces an extra pair of vector- like families 



[16 + 16 of SO(10)], at the TeV scale |3T], ||]. Adding such complete SO(10)-multiplets would 
of course preserve coupling unification. From the point of view of adding extra families, ESSM 
seems to be the minimal and also the maximal extension of the MSSM, that is allowed in that 
it is compatible with (a) LEP neutrino-counting, (b) precision electroweak tests, as well as (c) 
a semi-perturbative as opposed to non-perturbative gauge coupling unification |2l], |62| . The 



existence of two extra vector-like families of quarks and leptons can of course be tested at the LHC. 

Theoretical motivations for the case of ESSM arise on several grounds: (a) it provides a better 
chance for stabilizing the dilaton by having a semi-perturbative value for a un if ~ 0.35-0.3 [|2T|, in 
contrast to a very weak value of 0.04 for MSSM; (b) owing to increased two-loop effects |68|1 , 
it raises the unification scale Mx to (1/2 — 2) x 10 17 GeV and thereby considerably reduces the 
problem of a mismatch |28|] between the MSSM and the string unification scales (see Section 
[|); (c) It lowers the GUT-prediction for a 3 (mz) to (0.112-0.118) (in absence of unification-scale 
threshold corrections), which is in better agreement with the data than the corresponding value of 
(0.125-0.13) for MSSM; and (d) it provides a simple reason for inter-family mass-hierarchy |^T[ [d2 



In this sense, ESSM, though less economical than MSSM, offers some distinct advantages. 

In the present context, because of (b) and (c), ESSM naturally enhances the GUT-prediction 



for proton lifetime, in full accord with the data ||64|| . As explained in the appendix, the net 
result of these two effects — i.e. a raising of Mx and a lowering of ^(m^EssM — * s that for ESSM 
embedded in SO (10), tan (3 can span a wide range from 3 to even 30, and simultaneously the value 
or the upper limit on M e g- can range from (60 to 6) x 10 18 GeV, in full accord with our criterion 
for accurate coupling unification discussed above. 



19 Other authors (see e.g., Ref. |67]]) have considered proton decay in SUSY SO(10) by allowing for rather 
large GUT-scale threshold corrections, which do not, however, go well with our requirement of "natural coupling 
unification" . 

20 For instance, addition of two pairs of vector-like families at the TeV-scale, to the three chiral families, would 
cause gauge couplings to become non-perturbative below the unification scale. 
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As a result, in contrast to MSSM, ESSM allows for larger values of tan/3 (like 10 or 20), 
without needing large threshold corrections, and simultaneously without conflicting with the limit 
on proton lifetime. 

To be specific, consider first the case of a moderately large tan/3 = 10 (say), for which one 
obtains M e s « 1.8 x 10 19 GeV, with the "other" threshold correction — 5' 3 being about 5% (see 
Appendix for definition). In this case, one obtains: 



r-\i7K + ) std « 



(1.6-0.7) 
(10-1) 



{64 - 1/64} (7X10- years) • (49) 



As before, the upper and lower entries correspond to cases I (e' = 2 x 10~ 4 ) and II (e' = 0) of the 
fermion mass-matrix [see Eq. (|34])]. The uncertainty in the upper and lower entries in the square 
bracket of Eq. (|49p corresponds to that in the relative phases of the different contributions for the 
cases I and II respectively, while the factor {64-1/64} corresponds to uncertainties in the SUSY 
spectrum and the matrix element (see Eq. (|ID"D). 

We see that by allowing for an uncertainty of a factor of (30 — 100) jointly from the two 
brackets proton lifetime arising from the standard operators would be expected to lie in the range 
of (2.1 — 7) x 10 33 years, for the case of ESSM embedded in SO(10), even for a moderately large 
tan/3 = 10. Such a range is compatible with present limits, but accessible to searches in the near 
future. 

The other most important feature of ESSM is that, by allowing for larger values of M e g, 
especially for smaller values of tan/3 ~ 3 to 5 (say), the contribution of the standard operators by 
itself can be perfectly consistent with present limit on proton lifetime even for almost central or 
"median" values of the parameters pertaining to the SUSY spectrum, the relevant matrix element, 
e' and the phase- dependent factor. 

For instance, for ESSM, one obtains M eff « (4.5 x 10 19 GeV)(4/tan/3), with the "other" 
threshold correction - 5' 3 being about 5% [see Appendix and Eq. QAq) ]. Now, combining cases I 
(e' = 2 x 10 -4 ) and II (e' = 0), we see that the square bracket in Eq. ( f49|) which we will denote by 
[S], varies from 0.7 to 10, depending upon the relative phases of the different contributions and the 
values of e' . Thus as a "median" value, we will take [Sjmed ~ 2 to 6. The curly bracket {64-1/64}, 
to be denoted by {C}, represents the uncertainty in the SUSY spectrum and the matrix element 
[see Eq. ([40"D1. Again as a "nearly central" or "median" value, we will take {C} me d ~ 1/6 to 6. 
Setting M eff as above we obtain 

r- 1 (z7^ + ); t 7 dian " « [S] mcd {C} mod (0.45 x 10 33 years)(4/tan/3) 2 (SO(10)/ESSM). (50) 

Choosing a few sample values of the effective parameters [S] and {C}, with low values of tan (3 = 3 
to 5, the corresponding values of T~ 1 (uK + ), following from Eq. fl50|) , are listed below in Table |l[ 
Note that ignoring contributions from the new d = 5 operators for a moment]^], the entries in 
Table |l] represent a very plausible range of values for the proton lifetime, for the case of ESSM 
embedded in SO(10), with tan/3 ~ 3 to 5 (say), rather than upper limits for the same. This is 
because they are obtained for "nearly central" or "median" values of the parameters represented 
by the values of [S] and {C}, as discussed above. For instance, consider the cases {C}=1 and 
{C}=l/2 respectively, which (as may be inferred from the table) can quite plausibly yield proton 

21 As I will discuss in the next section, we of course expect the new d = 5 operators to be important and 
significantly influence proton lifetime (see e.g. Table ^). Entries in Table [l] could still represent the actual expected 



values of proton lifetimes, however, if the parameter K defined in |6.l| (also see p.5[ ) happens to be unexpectedly 
small (< 1). 
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lifetimes in the range of (2 to 5)xl0 33 years Now {C}=1 corresponds, e.g., to (3h = 0.014 GeV 3 
(the central value of Ref. J54J) rriq = 1.2 TeV and m^/rriq = 1/6 [see Eq. (|40|) 1, while that 
of {C}=l/2 would correspond, for example, to (3h = 0.014 GeV 3 , with nig ~ 710 GeV and 
m^r/rriq « 1/6. In short, for the case of ESSM, with low values of tan (3 ~ 5 to 5 fsm/j, squark 
masses can be well below 1 TeV, without conflicting with present limit on proton lifetime. This 
feature is not permissible within MSSM embedded in SO (10). 



Table 1: Proton lifetime, based on contributions from only the standard 
operators for the case of ESSM embedded in SO (10), with parameters 
being in the "median" range. 



tan (3 = 3 
[S]=3 
{C}=l/2 to 4 


tan (3 = 3 
[S]=6 
{C}=l/2 to 1 


tan/? = 5 
[S]=5.4 
{C}=1 to 6 


tan (3 = 5 

[S]=6 
{C}=1 to 4 


r-\uK+)i d SSM « 

(1.2 to 10) x 10 33 yrs 


r- 1 (pir+)gg w , « 

(2.5 to 5) x 10 33 yrs 


r-\vK+)i d SSM « 

(1.6 to 10) x 10 33 yrs 


(1.8 to 7.3) x 10 33 yrs 



Thus, confining for a moment to the standard operators only, if ESSM represents low-energy 
physics, and if tan (3 is rather small (3 to 5, say), we do not have to stretch the uncertainties 
in the SUSY spectrum and the matrix elements to their extreme values (in contrast to the case 
of MSSM) in order to understand why proton decay has not been seen as yet, and still can be 



optimistic that it ought to be discovered in the near future, with a lifetime < 10 years. The 
results for a wider variation of the parameters are listed in Table |2|, where contributions of the 
new d = 5 operators are also shown. 

It should also be remarked that if in the unlikely event, all the parameters (i.e. (3h, (m 



/m, 



m, 



and the phase-dependent factor) happen to be closer to their extreme values so as to extend 
proton lifetime, and if tan f3 is small (~ 3 to 5, say) and at the same time the value of M e g is close 
to its allowed upper limit (see Appendix), the standard d = 5 operators by themselves would tend 
to yield proton lifetimes exceeding even (0.8 to 2.5) xlO 34 years for the case of ESSM, (see Eq. 
(|49| ) and Table |2|). In this case (with the parameters having nearly extreme values), however, as 
I will discuss shortly, the contribution of the new d = 5 operators related to neutrino masses [see 
Eq. (|36"D1, are likely to dominate and quite naturally yield lifetimes bounded above in the range 
of (1 — 10) x 10 33 years (see Section 6J5 and Table §). Thus in the presence of the new operators, 
the range of (10 33 — 10 34 ) years for proton lifetime is not only very plausible but it also provides 
a reasonable upper limit, for the case of ESSM embedded in SO(10). 



6.5 Contribution from the new d=5 operators 



As mentioned in Section |6.1| , for supersymmetric G(224)/SO(10), there very likely exists a new 
set of d = 5 operators, related to neutrino masses, which can induce proton decay [see Eq. 
fl4"2])]. The decay amplitude for these operators for the leading mode (which in this case is v i ^K + ) 
becomes proportional to the quantity P = {( f 33) u (16 h) /M}h 33 K/ (Mi 6 tan 7), where (^33)^ and 
h>33 are the effective couplings defined in Eqs. (|7|) and fllS]) respectively, and M 16 and tan 7 are 
defined in the Appendix. The factor K, defined by (^33)^ = (f 3 3)i,K, is expected to be of order 
unity (see Section |0] for the origin of K). As a plausible range, we would take K 1/5 to 
2. Using M 16 tan7 = A'(16#) (see Appendix), and h 33 rs 1/2 (given by top mass), one gets: 
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P « [(f 33 ) u /M](l/2X')K. Here M denotes the string or the Planck scale (see Section [| and 
footnote 2); thus M « (1/2 — 1) X 10 18 GeV; and A' is a quartic coupling defined in the appendix. 
Validity of perturbative calculation suggests that A' should not much exceed unity, while other 



considerations suggest that A' should not be much less than unity either (see Ref. |14| , Section 6 
E). Thus, a plausible range for A' is given by A' ~ (1/2 — y/2). (Note it is only the upper limit 
on A' that is relevant to obtaining an upper limit on proton lifetime) . Finally, from consideration 
of v T mass, we have {f?s) v ~ 1 (see Section ||). We thus obtain: Ps;(5x 10 _19 GeV _1 )(l/\/2 to 
4)K. Incorporating a further uncertainty by a factor of (1/2 to 2) that arises due to choice of the 
relative phases of the different contributions (see Ref. [Q), the effective amplitude for the new 
operator is given by 

A(u II K + ) ncw « (1.5 x 10- 24 GeV" 1 )(l/2 v / 2 to 8)K (51) 

Note that this new contribution is independent of M e ff ; thus it is the same for ESSM as it is for 
MSSM, and it is independent of tan (3. Furthermore, it turns out that the new contribution is 
also insensitive to e'; thus it is nearly the same for cases I and II of the fermion mass-matrix. 
Comparing Eq. ([51]) with Eq. (f46j) we see that the new and the standard operators are typically 
quite comparable to one another. Since there is no reason to expect near cancellation between 
them (especially for both v^K^ and u T K + modes), we expect the net amplitude (standard+new) 
to be in the range exhibited by either one. It is thus useful to obtain the inverse decay rate 
assuming as if the new operator dominates. Substituting Eq. (|5l|) into Eq. (|4l|) and allowing for 
the presence of the v T K + mode with an estimated branching ratio of nearly 0.4 (see Ref. 
one obtains 

r^ 1 (77A' + ) llcw ss (0.25 x 10 31 years) [8-1/64] {64- 1/64} (K' 2 « 25 to 1/4) . (52) 

The square bracket represents the uncertainty reflected in Eq. fl51"|), while the curly bracket 
corresponds to that in the SUSY spectrum and matrix element (Eq. (|40|)). Allowing for the net 
uncertainty factor at the upper end, arising jointly from the three brackets in Eq. ( |52D to be 1000 
to 4000 (say), which can be realized for plausible range of values of the parameters (see below), 
the new operators related to neutrino masses, by themselves, lead to a proton decay lifetime given 
by: 

T- 1 (uK + )lZ eT « (2.5- 10) x 10 33 years (SO(10) or string G(224)) (Indep. of tan/3) . (53) 

The superscript "upper" corresponds to estimated lifetimes near the upper end. For instance, 
taking the curly bracket in Eq. ( |52"D to be rs 8 to 16 (say) [corresponding for example, to (3h = 
0.010 GeV 3 , (m^/mq) m 1/12 and m$ ps (1 to 1.4)(1.2 TeV)], instead of its extreme value of 64, 
and setting the square bracket in Eq. (|52|) to be ~ 6, and K~ 2 « 20, which are quite plausible, 
we obtain: T~ 1 (b / K + ) ncw « (2.5 — 5) x 10 33 years; independently of tan/3, for both MSSM and 
ESSM. Proton lifetime for other choices of parameters, which lead to similar conclusion, are listed 
in Table |. 

It should be stressed that the standard d — 5 operators [mediated by the color-triplets in the 
10h of SO(10)] may naturally be absent for a string-derived G(224)-model (see e.g. Ref. |3(J and 



PjJ), but the new d = 5 operators, related to the Majorana masses of the RH neutrinos and the 
CKM mixings, should very likely be present for such a model, as much as for SO(10). These would 
induce proton decay Q. Thus our expectations for the proton decay lifetime [as shown in Eq. ( "55) J 

22 In addition, quantum gravity induced d = 5 operators are also expected to be present at some level, depending 
upon the degree of suppression of these operators due to flavor symmetries (see e.g. Ref. J34|). 
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and the prominence of the n + K° mode (see below) hold for a string- derived G (224) -model, i us ^ 
as they do for SO(10). For a string - G(224) - model, however, the new d=5 operators would be 
essentially the sole source of proton decay 21 . 

Nearly the same situation emerges for the case of ESSM embedded in G(224) or SO(10), with 
low tan/3(^ 3 to 10, say), especially if the parameters (including /3 H , m^/rrig, nig, the phase- 
dependent factor as well as M e g) happen to be somewhat closer to their extreme values so as 
to extend proton lifetime. In this case, (that is for ESSM) as noted in the previous sub-section, 
the contribution of the standard d = 5 operators would be suppressed; and proton decay would 
proceed primarily via the new operators with a lifetime quite plausibly in the range of 10 33 — 10 34 
years, as exhibited above. 



6.6 The Charged Lepton Decay Modes (p — ► p, + K° and p — > e + ir°) 

I now note a distinguishing feature of the SO(10) or the G(224) model presented here. Allowing 
for uncertainties in the way the standard and the new operators can combine with each other for 
the three leading modes i.e. T7 T K + , V^K + and fi + K°, we obtain (see Ref. []nj for details): 

B(fi + K°) std+new « [1% to 50%] k (SO(10) or string G(224)) (54) 

where k denotes the ratio of the squares of relevant matrix elements for the [i + K° and VK + modes. 
In the absence of a reliable lattice calculation for the uK + mode, one should remain open to the 
possibility of k « 1/2 to 1 (say). We find that for a large range of parameters, the branching ratio 
B(fx + K°) can lie in the range of 20 to 40% (if k ~ 1). This prominence of the fi + K° mode for 
the SO(10)/G(224) model is primarily due to contributions from the new d = 5 operators. This 
contrasts sharply with the minimal SU(5) model, in which the fi + K° mode is expected to have a 
branching ratio of only about 10 -3 . In short, prominence of the fi + K° mode, if seen, would clearly 
show the relevance of the new operators, and thereby reveal the proposed link between neutrino 
masses and proton decay . 

The d = 5 operators as described here (standard and new) would lead to highly suppressed 
e + 7r° mode, for MSSM or ESSM embedded in SO(10). The gauge boson-mediated d — 6 operators, 
however, still give (using the recently determined matrix element an = 0.015 ± 0.001 GeV 3 |54|) 
proton decaying into e + 7r° with an inverse rate: 



r- 1 ^ - e+^sT^ « 10 35±1 years . (55) 

This can well be as short as about 10 34 years. For the case of ESSM embedded into SO(10) [or 
for an analogous case embedded into SU(5)], there are two new features. Considering that in this 
case, both a un jf and the unification scale Mx (thereby the mass My of the (X, Y) gauge bosons) 
are raised by nearly a factor of (6 to 7) and (2.5 to 5) respectively, compared to those for MSSM 
(see discussions in Section p.4[ ), and that the inverse decay rate is proportional to (My/a^ nif ), we 
expect 

r-'d- - e V)S8r um - (i t° i7)r-'(P - ^ ) S M O S ,/su(5, ■ (56) 

The net upshot is that the gauge boson-mediated d = 6 operators can quite plausibly lead to 
observable e + ir° decay mode with an inverse decay rate in the range of 10 34 -10 35 years. For ESSM 
embedded in SO (10), there can be the interesting situation that both vK + (arising from d — 5) 
and e + 7r° (arising from d = 6) may have comparable rates, with proton having a lifetime ~ (1/2-2) 
x 10 34 years. It should be stressed that the e + 7r°-mode is the common denominator of all GUT 
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models (SU(5), SO(10), etc.) which unify quarks and leptons and the three gauge forces. Its 
rate as mentioned above is determined essentially by the SUSY unification-scale, without the 
uncertainty of the SUSY-spectrum. I should also mention that the e + 7r°-mode is predicted to be 
the dominant mode in the flipped SU(5) x U(l)-model |){J. For these reasons, intensifying the 
search for the e + 7r°-mode to the level of sensitivity of about 10 35 years in the next generation 
proton decay detector should be well worth the effort. 

Before summarizing the results of this section, I note below a few distinctive features of the 
conventional approach adopted here compared to those of some alternatives. 



6.7 Conventional Versus Other Approaches 

In these lectures, as elaborated in Section 3, I have pursued systematically the consequences for 
fermion masses, neutrino oscillations and proton decay of the assumption that essentially the con- 
ventional picture of SUSY grand unification |3|, |], [5], [6|, [7] holds, providing a good effective theory 
in 4D between the conventional GUT-scale M x ~2x 10 16 GeV (for ESSM, M x ~ (1/2-2) x 10 17 
GeV) and the conventional string scale M st ~ (few to 10) x 10 17 GeV. Believing in an underly- 
ing string/M-theory, and yet knowing that a preferred ground state of this theory is not yet in 
hand, the attitude, based on a bottom-up approach, has been to subject the assumed effective 
theory of grand unification to as many low-energy tests as possible, and to assess its soundness 
on empirical grounds. With this in mind, I have assumed that either a realistic 4D SO (10)- 
solution (with the desired mechanism of doublet-triplet splitting operating in 4D), or a suitable 
string-derived G(224)-solution (with Mx ~ (l/2)M st , see footnote 2) emerges effectively from 
an underlying string theory at the conventional string scale as mentioned above, and that the 
G(224)/SO(10) symmetry breaks into G(213) at the conventional GUT-scale Mx- The extra 
dimensions of string/M-theory are assumed to be tiny lying between the GUT-scale size ~ M^ 1 
and the string-size M^ , so as not to disturb the successes of GUT (see below). As mentioned 
before, this conventional picture of grand unification described above seems to be directly moti- 
vated on observational grounds such as those based on (a) coupling unification or equivalently the 
agreement between the observed and the predicted values of sin 2 9w (see Section 3), (b) neutrino 
masses including Am 2 (z/ /J -i/ T ) and (c) the fact that spontaneous violation of B L local symmetry 
seems to be needed to implement baryogenesis via leptogenesis [ P~5] , [T7| . The relevance of the 



group theory of G(224)/SO(10)-symmetry for the 4D theory is further suggested by the success 
of the predictions of the masses and the mixings of all fermions including neutrinos; these include 
m" ~ m°, ^(^Dirac) ~ m t(Mx), and the smallness of V c b ~ 0.04 correlated with the largeness of 
sin 2 20°^ « 1 (see Section 5). 

In contrast to this conventional approach based on a presumed string- unified G(224) or an 
SO(10)-symmetry, there are several alternative approaches (scenarios) which have been proposed 
in the literature in recent years. Of importance is the fact that in many of these alternatives an 
attempt is made to strongly suppress proton decay, in some cases exclusively the d = 5 operators 
(though not necessarily the d = 6), invariably utilizing a higher dimensional mechanism. Each 
of these alternatives is interesting in its own right. However, it seems to me that the collection 
of successes mentioned above is not (yet) realized within these alternatives. For comparison, I 
mention briefly only a few, leaving out many interesting variants. 

One class of alternatives is based on the idea of large extra dimensions |36| , [38|j , or a low 
string scale of order one TeV |[7). Though most intriguing, it does not seem to provide simple 



explanations for (a) coupling unification, (b) neutrino- masses (or their (mass) 2 -differences) of the 
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observed magnitudes^], (c) a large (or maximal) v^-Vt oscillation angle, and (d) baryogenesis via 
leptogenesis that seems to require violation of B L at high temperatures. Within this scenario, 
quantum-gravity induced proton decay would ordinarily be extra rapid. This is prevented, for 
example, by assuming that quarks and leptons live in different positions in the extra dimension. 
It appears to me that this idea (introduced just to prevent proton decay) however, sacrifices the 
simple reason for the co-existence quarks and leptons that is provided by a gauge unification of 
matter within a family as in G(224) or SO(10). 

There is an alternative class of attempts, carried out again in the context of higher dimensional 
theories, which, in contrast to the case mentioned above, assume that the extra dimensions (d > 4) 
are all small, lying between (or around) the conventional GUT and string scales. The approach 
of this class of attempts is rather close in spirit to that of the conventional approach of grand 
unification pursued here (see Section 3). As may be seen from the discussions below, they could 
essentially coincide with the string- unified G(224)-picture presented here if the effective symmetry 
in 4D, below the string (or compactification) scale, contains at least the G(224) symmetry. 

Motivated by the original attempts carried out in the context of string theory |n] most of the 
recent attempts in the class mentioned above are made in the spirit of a bottom-up approach^] 
to physics near the GUT and the string scales. They assume, following the spirit of the results 
of Ref. ||71|| , and of analogous results obtained for the free fermionic formulation of string theory 



|72f (for applications based on this formulation, see e.g., Ref. f73fl , p0| , |3~Ifl and [[7§]), that 
grand unification occurs, through symmetries like Eg, SO(10) or SU(5), only in some higher 
dimension (d > 4), and that the breaking of the unification gauge symmetry to some lower 
symmetry containing the standard model gauge group as well as doublet-triplet splitting occurs 
in the process of compactification. More specifically the latter two phenomena take place through 
either (a) Wilson lines f7lfl , or (b) orbifolds |75| (for an incomplete list of recent attempts based 



on orbifold compactification, see e.g., Refs. [76], \n\ |78], |79], |80], |82|, |3|, |J, pq|), or (c) essentially 
equivalently by a set of boundary conditions together with the associated GSO projections for 
the free fermionic formulation (see e.g., f3(J, |31, [73], 74]]), or (d) discrete symmetries operating in 



higher dimensions [56 



Most of these attempts end up not only in achieving (a) doublet-triplet splitting by projecting 
out the relevant color triplets from the zero mode-spectrum in 4D, and (b) gauge symmetry 
breaking, as mentioned above, but also (c) suppressing strongly or eliminating the d = 5 proton 
decay operators. It should be mentioned, however, that in some of these attempts (see e.g., |77|). 



the mass of the X gauge boson is suggested to be lower than the conventional GUT-scale of 
2xl0 16 GeV by about a factor of 3 to 8; correspondingly they raise the prospect for observing the 
d = 6 gauge boson mediated e + n° mode, which is allowed in [[F7[ . 

One crucial distinction between the various cases is provided by the nature of the effective 
gauge symmetry that is realized in 4D, below the string (or compactification) scale. References 



[ Tq , |77| , |7q , |80| , pi] , p2|] assume a supersymmetric SU(5) gauge symmetry in 5D, which is bro- 



ken down to the standard model gauge symmetry in 4D through compactification. References 



83 1 and |84| , on the other hand, assume a supersymmetric SO(10) gauge symmetry in 6D and 
show (interestingly enough) that there are two 5D subspaces containing G(224) and SU(5)xU(l) 
subgroups respectively, whose intersection leads to SU(3)xSU(2)xU(l)yxU(l)x in 4D, which 



23 By placing the singlet (right-handed) neutrino in the bulk, for example, one can get a light Dirac neutrino |70|| 
with a mass m„ w kv ew M*/M pi w k(2 x 10~ 5 eV), where M* w 1 TeV, M Pl w 10 19 GeV (as in JTo) ) , and k is the 



effective Yukawa coupling. To get m v ~ 1/20 eV (for SuperK), one would, however, need too large ac~2x 10 
and/or too large a value for M* > 100 TeV; which would seem to face the gauge-hierarchy problem. 
24 This is of course also the case for the approach adopted here which is outlined in Section 3. 
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contains B-L. [An alternative construction using SO(10) in 5D leads just to the standard model 



gauge group in 4D [85].] While it is desirable to have B-L in 4D, consistent breaking of U(l)x 
(or B-L) and generating desired masses of the right handed neutrinos, not to mention the masses 
and the mixings of the other fermions, is not yet realized in these constructions. 

For comparison, it seems to me that at the very least B-L should emerge as a generator in 4D to 
implement baryogenesis via leptogenesis, and also to protect RH neutrinos from acquiring a string- 
scale mass. This feature is not available in models which start with SU(5) in 5D. Furthermore, 
the full SU(4)-color symmetry, which of course contains B-L, plays a crucial role in yielding not 
only m' rs mj but also (a) ^(^Dirac) ~ m t(Mx) that is needed to account for m{y T ) or rather 
Am 2 (z/^-i/ T ), in accord with observation (see Section 4), and (b) the smallness of V c b together with 
the near maximality of sin 2 29° sc v _ r (see Section 5). The symmetry SU(2) L xSU(2) R is also most 
useful in that it relates the masses and mixings of the up and the down sectors. Without such 
relations, we will not have the predictivity of the framework presented in Section 5. 

In short, as mentioned before, certain intriguing features of the masses and mixings of all 
fermions including neutrinos, of the type mentioned above, as well as the need for leptogenesis, 
seem to strongly suggest that the effective symmetry below the string-scale in 4D should con- 
tain minimally the symmetry G(224) [or a close relative G(214)] and maximally SO(10). The 
G(224)/SO(10)-framework developed here has turned out to be the most predictive, in large part 
by virtue of its group structure and the assumption of minimality of the Higgs system. Given 
that it is also most successful so far, as regards its predictions, derivation of such a picture from 
an underlying theory, especially at least that based on an effective G (224)-symmetryp s | in 4D 
leading to the pattern of Yukawa couplings presented here remains a challenge.^ Pending such a 
derivation, however, given the empirical support it has received so far, it makes sense to test the 
supersymmetric G(224)/SO(10)-framework, and thereby the conventional picture of grand unifi- 
cation on which it rests, thoroughly. There are two notable missing pieces of this picture. One is 
supersymmetry which will be probed at the LHC and a future NLC. The other, that constitutes 
the hallmark of grand unification, is proton decay. The results of this section on proton decay are 
summarized below. 



6.8 Section Summary 

Given the importance of proton decay, a systematic study of this process has been carried out 
within the supersymmetric SO(10)/G(224)-frameworkFJ, with special attention paid to its depen- 
dence on fermion masses and threshold effects. A representative set of results corresponding to 
different choices of parameters is presented in Tables [TJ and Allowing for the ESSM- variant, the 

25 For this case, following the examples of Refs. |3(]] and fl3l| , the color triplets in the 10/r of SO(10) would be 
projected out of the zero-mode spectrum, and thus the standard d — 5 operators which would have been induced 
by the exchange of such triplets would be absent, as in Refs. |7(| [77| (78[ |8(J |8l], ||[ |||, |84|, ^6). But, as long 
as the Majorana masses of the RH neutrinos are generated as in Section 4, the new neutrino-mass related d = 5 
proton decay operators would generically be present (see Section 6 E). 

26 In this regard, three-generation solutions containing the G(224)-symmetry in 4D have been obtained in the 
context of the fcrmionic formulation of string theory in Ref. ]30[| , within type-I string vacua with or without 
supersymmetry in [^7[ |8[ in the context of D-brane inspired models in |9(J , within type-I string-construction 
or string-motivated models obtained from intersecting D-branes (with G(224) breaking into G(213) at Mx ~ M st ) 
in J9l], ^2), in string model with unification at the string scale in |93|, and in other contexts (see e.g. [Q and [95]). 

27 As described in Sections ||, |] and |[ 
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study strongly suggests that an upper limit on proton lifetime is given by 

Voton < (1/3-2) x 10 34 years, (57) 

with VK + being the dominant decay mode, and quite possibly fi + K° and e + 7r° being prominent. 
Although there are uncertainties in the matrix element, in the SUSY-spectrum, in the phase- 
dependent factor, tan {3 and in certain sensitive elements of the fermion mass matrix, notably e' 
(see Eq. fl48|) for predictions in cases I versus II), this upper limit is obtained, for the case of 
MSSM embedded in SO(10), by allowing for a generous range in these parameters and stretching 
all of them in the same direction so as to extend proton lifetime. In this sense, while the predicted 
lifetime spans a wide range, the upper limit quoted above, in fact more like 10 33 years, is most 
conservative, for the case of MSSM (see Eq. (|48| ) and Table [[]). It is thus tightly constrained 
already by the empirical lower limit on r~ 1 (i7.K' + ) of 1.9 x 10 33 years to the point of being disfa- 
vored. For the case of ESSM embedded in SO(10), the standard d — 5 operators are suppressed 
compared to the case of MSSM; as a result, by themselves they can naturally lead to lifetimes in 
the range of (1 — 10) x 10 33 years, for nearly central values of the parameters pertaining to the 
SUSY-spectrum and the matrix element (see Eq. (|50|) and Table ||). Including the contribution 
of the new d = 5 operators, and allowing for a wide variation of the parameters mentioned above, 
one finds that the range of (10 33 — 2 x 10 34 ) years for proton lifetime is not only very plausible 
but it also provides a rather conservative upper limit, for the case of ESSM embedded in either 
SO(10) or G(224) (see Section [67| and Table g). Thus our study provides a clear reason to expect 
that the discovery of proton decay should be imminent for the case of ESSM, and even more so 
for that of MSSM. The implication of this prediction for a next-generation detector is emphasized 
in the next section. 

7 Concluding Remarks 

The preceding sections show that, but for two missing pieces — supersymmetry and proton decay — 
the evidence in support of grand unification is now strong. It includes: (i) the observed family- 
structure, (ii) quantization of electric charge, (iii) the meeting of the gauge couplings, (iv) neutrino- 
oscillations as observed at SuperK, (v) the intricate pattern of the masses and mixings of all 
fermions, including the neutrinos, and (vi) the need for B L as a generator, to implement baryo- 
genesis. Taken together, these not only favor grand unification but in fact select out a particular 
route to such unification, based on the ideas of supersymmetry, SU(4)-color and left-right sym- 
metry. Thus they point to the relevance of an effective string-unified G(224) or SO(10)-symmetry 
in four dimensions, as discussed in Sections 3 and 4. 

Based on a systematic study of proton decay within the supersymmetric SO(10)/G(224)- 
framework, that (a) allows for the possibilities of both MSSM and ESSM, and (b) incorporates 
the improved values of the matrix element and renormalization effects, I have argued that a 
conservative upper limit on the proton lifetime is about (1/3-2) xlO 34 years. 

So, unless the fitting of all the pieces listed above is a mere coincidence, it is hard to believe 
that that is the case, discovery of proton decay should be around the corner. In particular, as 
mentioned in the Introduction, one expects that candidate events should very likely be observed 
in the near future already at SuperK, if its operation is restored. However, allowing for the 
possibility that proton lifetime may well be near the upper limit stated above, a next-generation 
detector providing a net gain in sensitivity by a factor five to ten, compared to SuperK, would be 
needed to produce real events and distinguish them unambiguously from the background. Such 
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an improved detector would of course be essential to study the branching ratios of certain crucial 
though (possibly) sub-dominant decay modes such as the fj, + K° and e + n° as mentioned in Section 

The reason for pleading for such improved searches is that proton decay would provide us 
with a wealth of knowledge about physics at truly short distances (< 10 -30 cm), which cannot 
be gained by any other means. Specifically, the observation of proton decay, at a rate suggested 
above, with T>K + mode being dominant, would not only reveal the underlying unity of quarks 
and leptons but also the relevance of supersymmetry. It would also confirm a unification of the 
fundamental forces at a scale of order 2 x 10 16 GeV. Furthermore, prominence of the fi + K° mode, 
if seen, would have even deeper significance, in that in addition to supporting the three features 
mentioned above, it would also reveal the link between neutrino masses and proton decay, as 
discussed in Section |J In this sense, the role of proton decay in probing into physics at the most 
fundamental level is unique. In view of how valuable such a probe would be and the fact that the 
predicted upper limit on the proton lifetime is at most a factor of three to ten higher than the 
empirical lower limit, the argument in favor of building an improved detector seems compelling. 

To conclude, the discovery of proton decay would undoubtedly constitute a landmark in the 
history of physics. It would provide the last, missing piece of gauge unification and would shed 
light on how such a unification may be extended to include gravity in the context of a deeper 
theory. 
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APPENDIX: A Natural Doublet-Triplet Splitting Mechanism in SO (10) 

In supersymmetric SO(10), a natural doublet-triplet splitting can be achieved by coupling the 
adjoint Higgs 45h to a 10h and a 10^, with 45h acquiring a unification-scale VEV in the 
B-L direction |96|, j97fl : (45h) = (a, a, a, 0, 0) x r 2 with a ~ My. As discussed in Section ^ 
to generate CKM mixing for fermions we require (16h)<2 to acquire a VEV of the electroweak 
scale. To ensure accurate gauge coupling unification, the effective low energy theory should not 
contain split multiplets beyond those of MSSM. Thus the MSSM Higgs doublets must be linear 
combinations of the SU(2)x doublets in 10h and 16h- A simple set of superpotential terms that 
ensures this and incorporates doublet-triplet splitting is [ [Hfl : 

W H = A10h45h10' h + M m 10' u 2 + A'TetfTez/lOH + M W 16 H 16 H . (Al) 

A complete superpotential for 45h, 16h, 16h, 10#, 10^ and possibly other fields, which ensure 
that (a) 45h, 16h and 16h acquire unification scale VEVs with (45h) being along the (B-L) 
direction; (b) that exactly two Higgs doublets (H u , H^) remain light, with being a linear 
combination of (IOh)^ and (16n)d', and (c) there are no unwanted pseudoGoldstone bosons, can 
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be constructed. With (45h) in the B L direction, it does not contribute to the Higgs doublet 
mass matrix, so one pair of Higgs doublet remains light, while all triplets acquire unification scale 
masses. The light MSSM Higgs doublets are [13] 



H u = 10 u , H d = cos 7 10 d + sin 7 16 d , (A2) 

with tan 7 = A'(16 H )/M 16 . Consequently, (10} d = (cos~/)v d , (lQ d ) = (siwy)v d , with (H d ) = v d 
and (16^) and (10,^) denoting the electroweak VEVs of those multiplets. Note that H u is purely 
in 10h and that (lO^) 2 + (16 d ) 2 = v\. This mechanism of doublet-triplet (DT) splitting is the 
simplest for the minimal Higgs systems. It has the advantage that it meets the requirements of 
both D-T splitting and CKM-mixing. In turn, it has three special consequences: 

(i) It modifies the familiar SO(10)-relation tan/3 = v u /v d = nit/mf, « 60 to Q: 

tan/3/ cos7 ~ m t /mi, ~ 60 . (A3) 

As a result, even low to moderate values of tan/5 ss 3 to 10 (say) are perfectly allowed in SO(10) 
(corresponding to cos 7 ss 1/20 to 1/6). 

(ii) The most important consequence of the DT-splitting mechanism outlined above is this: 
In contrast to SU(5), for which the strengths of the standard d = 5 operators are proportional 
to (My,) -1 (where M Hc ~ few x 10 16 GeV (see Eq. (!§), for the SO(10)-model, they become 
proportional to M~g, where M cS = (\a) 2 /M 10 > ~ M x /M 10 >. As noted in Ref. [0, M io> can 



be naturally smaller (due to flavor symmetries) than Mx and thus M e g correspondingly larger 
than Mx by even one to three orders of magnitude. Now the proton decay amplitudes for SO(10) 
in fact possess an intrinsic enhancement compared to those for SU(5), owing primarily due to 
differences in their Yukawa couplings for the up sector (see Appendix C in Ref. |L4|]). As a result, 
these larger values of M eff ~ (10 18 — 10 19 ) GeV are in fact needed for the SO(10)-model to be 
compatible with the observed limit on the proton lifetime. At the same time, being bounded 
above by considerations of threshold effects (see below), they allow optimism as regards future 
observation of proton decay. 

(iii) M e fj gets bounded above by considerations of coupling unification and GUT-scale threshold 
effects as follows. Let us recall that in the absence of unification-scale threshold and Planck-scale 
effects, the MSSM value of a 3 {m z ) in the MS scheme, obtained by assuming gauge coupling 
unification, is given by «3(^z)mssm = 0.125 — 0.13 [0. This is about 5 to 8% higher than 



the observed value: a 3 (m z ) = 0.118 ± 0.003 []13| . Now, assuming coupling unification, the net 



(observed) value of a 3} for the case of MSSM embedded in SU(5) or SO(10), is given by: 

a 3 (m z ) nct = a 3 (m z )° MSSM + Aa 3 (m z )^ SM + A' 3 (A4) 

where Aa 3 (m^)DT and A' 3 represent GUT-scale threshold corrections respectively due to doublet- 
triplet splitting and the splittings in the other multiplets (like the gauge and the Higgs multiplets), 
all of which are evaluated at m z . Now, owing to mixing between 10^ and lQ d [see Eq. (|A2|) 1, one 
finds that Aa 3 (m z ) BT is given by [a 3 (m z ) 2 /2ir](9/7) ln(M efr cos 7 /Mx) 0. 



As mentioned above, constraint from proton lifetime sets a lower limit on M e fj given by M e g > 
[1 — 6) x 10 18 GeV. Thus, even for small tan/3 w 2 (i.e. cos 7 w tan(/3/60) w 1/30), Aa 3 (m z ) DT i s 



28 It is worth noting that the simple relationship between cos 7 and tan/3 — i.e. cos 7 « tan (3/(mt/mb) — would 
be modified if the superpotential contains an additional term like A"16# • 16# • 10^, which would induce a mixing 
between the doublets in 10^, 16^ and 10^. That in turn will mean that the upper limit on M c ff COS7 following 
from considerations of threshold corrections (see below) will not be strictly proportional to tan/3. I thank Kaladi 
Babu for making this observation. 



36 



positive; and it increases logarithmically with M e ff. Since «3(^z)mssm i s higher than a 3 (mz)obs, 
and as we saw, Aa 3 (mz)-DT is positive, it follows that the corrections due to other multiplets 
denoted by 6' 3 = A 3 /a 3 (mz) should be appropriately negative so that a 3 {mz)net would agree with 
the observed value. 

In order that coupling unification may be regarded as a natural prediction of SUSY unification, 
as opposed to being a mere coincidence, it is important that the magnitude of the net other 
threshold corrections, denoted by 5' 3 , be negative but not any more than about 8 to 10% in 
magnitude (i.e. —5' 3 < (8 — 10)%). It was shown in Ref. [14] that the contributions from the 



gauge and the minimal set of Higgs multiplets (i.e. 45#, 16#, 16# and 10#) leads to threshold 
correction, denoted by 5' 3 , which has in fact a negative sign and quite naturally a magnitude of 4 
to 8%, as needed to account for the observed coupling unification. The correction to a 3 {mz) due 
to Planck scale physics through the effective operator F^F^A^h/M does not alter the estimate 
of 5' 3 because it vanishes due to antisymmetry in the SO(10)- contraction. 

Imposing that 5' 3 (evaluated at m^)be negative and not any more than about 10-11% in 
magnitude in turn provides a restriction on how big the correction due to doublet-triplet splitting— 
i.e. Aa 3 (mz)Df — can be. That in turn sets an upper limit on M eff cos7, and thereby on M eff 
for a given tan/5. For instance, for MSSM, with tan/3 = (2,3,8), one obtains (see Ref. [|14| ): 
M e ff < (4,2.66, 1) x 10 18 GeV. Thus, conservatively, taking tan/3 > 3, one obtains: 

M eS < 2.7 x 10 18 GeV (MSSM) (tan/3>3). (A5) 



Limit on M eff For The case of ESSM 

Next consider the restriction on M e g that would arise for the case of the extended supersymmetric 
standard model (ESSM), which introduces an extra pair of vector-like families (16 + 16) of SO(10)) 
at the TeV scale [^](see also footnote 16). In this case, a un ff is raised to 0.25 to 0.3, compared 
to 0.04 in MSSM. Owing to increased two-loop effects the scale of unification Mx is raised to 
(1/2 - 2) x 10 17 GeV, while a 3 (m z )^ SSM is lowered to about 0.112-0.118 pT], |8|. 

With raised Mx, the product M e ff cos7 ~ M e ff(tan/3)/60 can be higher by almost a factor of 
five compared to that for MSSM, without altering Aa 3 (mz)DT- Furthermore, since ce 3 (mz)MssM 
is typically lower than the observed value of a 3 {mz) (contrast this with the case of ESSM), for 
ESSM, M e ff can be higher than that for MSSM by as much as a factor of 2 to 3, without requiring 
an enhancement of 5' 3 . The net result is that for ESSM embedded in SO(10), tan/3 can span a 
wide range from 3 to even 30 (say) and simultaneously the upper limit on M e s can vary over the 
range (60 to 6)xl0 18 GeV^, satisfying 

M cfr < (6 x 10 18 GeV)(30/tan/3) (ESSM), (A6) 

with the unification-scale threshold corrections from "other" sources denoted by S 3 = A 3 /a 3 (mz) 
being negative, but no more than about 5% in magnitude. As noted above, such values of 5' 3 
emerge quite naturally for the minimal Higgs system. Thus, one important consequence of ESSM 
is that by allowing for larger values of M e g (compared to MSSM), without entailing larger values 
of 5' 3 , it can be perfectly compatible with the limit on proton lifetime for almost central values of 
the parameters pertaining to the SUSY spectrum and the relevant matrix elements (see Eq. (40)). 
Further, larger values of tan/? (10 to 30, say) can be compatible with proton lifetime only for the 
case of ESSM, but not for MSSM. These features are discussed in the text, and also exhibited in 
Table |. 
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Table 2: Values of proton lifetime (T" 
parameters. 



1 (p — > vK + ))for a wide range of 



Parameters 
(spectrum/Matrix 


MSSM - 
Std. 


+ SO(10) 
d=5 


ESSM - 
Std. 


► SO(10) 
d=5 


f MSSM 1 

{ or } -> G(224)/SO(10) 
[ ESSM J 


element) 


Intermed. 


s' & phase''' 


Intermed. e' & phase^ 


New d=5 ft 




tan/3=3 


tan/3=10 


tan (3=5 


tan p=10 


Independent of tan /? 


Nearly "central" 


0.2 x 10 32 


1.6xl0 30 


0.25xl0 34 


0.7xl0 33 


0.50 x 10 33 


{C}=2 


yrs 


yrs 


yrs* 


yrs 


yrs'"" 


Intermediate 


0.7xl0 32 


0.6xl0 31 


lxlO 34 


2.8xl0 33 


2xl0 33 


{C}=8 


yrs 


yrs 


yrs* 


yrs 


yrs'"' 


Nearly Extreme 


0.3 xlO 33 


2.6xl0 31 


4xl0 34 


l.lxlO 34 


8xl0 33 


{C}=32 


yrs 


yrs 


yrs* 


yrs 


yrs'"'" 



*In this case, lifetime is given by the last column. 

• Since we are interested in exhibiting expected proton lifetime near the upper end, we are 
not showing entries in Table H corresponding to values of the parameters for the SUSY spectrum 
and the matrix element [see Eq. 



i0|) , for which the curly bracket {C} appearing in Eqs. (47Q, 
would be less than one (see however Table |l|). In this context, we have chosen here 
"nearly central", "intermediate" and "nearly extreme" values of the parameters such that the 
said curly bracket is given by 2, 8 and 32 respectively, instead of its extreme upper-end value 



of 64. For instance, the curly bracket would be 2 if (3h 

(3 H = 0.010 
0.007 GeV 3 



m w /rriq « (1/7.2), while it would be 8 if H = 0.010 GeV 3 , 

"' ■"' ~ " ' m 9 - ~ 



(0.0117) GeV 3 , 



rriq ~ 



1.44 TeV and m-^/m 



1.2 TeV and 
,~1/10; 
1/12. 



and it would be 32 if, for example, (3h = 0.007 GeV 3 , rriq ~ \/2(L2 TeV) and m^/rriq 

f All the entries for the standard d = 5 operators correspond to taking an intermediate value 
of e' ~ (1 to 1.4) x 10~ 4 (as opposed to the extreme values of 2 x 10 -4 and zero for cases I and 
II, see Eq. (|34])) and an intermediate phase-dependent factor such that the uncertainty factor in 
the square bracket appearing in Eqs. (|7|) and ( |49|) is given by 5, instead of its extreme values of 
2x4 = 8 and 2.5 x 4 = 10, respectively. 

ft For the new operators, the factor [8-1/64] appearing in Eq. ( |52j) is taken to be 6, and K~ 2 , 
defined in Section p.l\ , is taken to be 25, which are quite plausible, in so far as we wish to obtain 
reasonable values for proton lifetime at the upper end. 

• The standard d = 5 operators for both MSSM and ESSM are evaluated by taking the 
upper limit on M eS (defined in the text) that is allowed by the requirement of natural coupling 
unification. This requirement restricts threshold corrections and thereby sets an upper limit on 
M eff , for a given tan/3 (see Section || and Appendix). 

• For all cases, the standard and the new d = 5 operators must be combined to obtain the net 
amplitude. For the three cases of ESSM marked with an asterisk, and other similar cases which 
arise for low tan (3 ~ 3 to 6 (say), the standard d = 5 operators by themselves would lead to proton 
lifetimes typically exceeding (0.25-4) x 10 34 years. For these cases, however, the contribution from 
the new d = 5 operators would dominate, which quite naturally lead to lifetimes in the range of 
(10 33 — 10 34 ) years (see last column). 

• As shown above, the case of MSSM embedded in SO(10) is tightly constrained to the point 
of being disfavored by present empirical lower limit on proton lifetime Eq. fl4"2"D [see discussion 
following Eq. (M)]. 
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• Including contributions from the standard and the new operators, the case of ESSM, embed- 
ded in either G(224) or SO(10), is, however, fully consistent with present limits on proton lifetime 
for a wide range of parameters; at the same time it provides optimism that proton decay will be 
discovered in the near future, with a lifetime < 10 34 years. 

• The lower limits on proton lifetime are not exhibited. In the presence of the new operators, 
these can typically be as low as about 10 29 years (even for the case of ESSM embedded in SO(10)). 
Such limits and even higher are of course long excluded by experiments. 
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